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—When a fluid flows through a vessel at a constant rate, either “ piston-flow " or perfect 
mixing is usually assumed. In practice many systems do not conform to either of these assumptions, 
so that calculations based on them may be inaccurate. It is explained how distribution-functions for 
residence-times can be defined and measured for actual systems. Open and packed tubes are discussed 
as systems about which predictions can be made. The use of the distribution-functions is illustrated by 
showing how they can be used to calculate the efficiencies of reactors and blenders. It is shown how 
models may be used to predict the distribution of residence-times in large systems. 


Résumé—Quand, dans un récipient, on introduit, 4 vitesse constante, un fluide donné, on suppose 
généralement soit un mélange parfait, soit un “* écoulement frontal parfait.” En pratique, de nombreux 
systémes s’écartent de l'une ou l'autre de ces hypothéses simplificatrices et les calculs qui en résultent 
sont plus ou moins inexacts. L’auteur expose, pour des syst¢mes réels, comment l'on peut définir et mesurer 
des fonctions de distribution pour la “ durée de séjour ™: ceci peut s‘appliquer 4 des tubes vides ou munis 
de garnissages. Par emploi de ces fonctions de distribution, l'auteur montre comment on peut calculer 
lefficacité des réacteurs ou des mélangeurs. Des modéles peuvent étre utilisés pour prévoir la répartition 


When a stream of material flows steadily through a 
vessel such as a pipe or a tank, in which it takes part 
in some process such as chemical reaction, heat- or 
mass-transfer, or simple mixing, it is usual to make 
use of one of the following assumptions for the pur- 
poses of calculation : 

(a) The fluid in the vessel is completely mixed, so 
that its properties are uniform and identical with 
those of the outgoing stream. This assumption is 
frequently made the basis of calculations on stirred 
reactors or blenders. 

(b) Elements of fluid which enter the vessel at the 
same moment move through it with constant and 
equal velocity on parallel paths, and leave at the 
saine moment. This type of behaviour will be referred 
to as “ piston flow,”’ and is normally assumed when 
considering flow through heat-exchangers, catalytic 
reactors, packed towers, chromatographic columns, 
etc. 
It is clear that there are many cases in which 
neither type of flow corresponds exactly to the facts - 
for instance, fluid in laminar flow in a pipe, or gas 


des “* durées de séjour dans des systémes de grandes dimensions. 


flowing through a fluidised catalytic reactor or blast 
furnace. It is of some importance to investigate the 
discrepancies between the assumed and actual 
behaviour of such systems, and where necessary to 
allow for them in making calculations. This paper 
is intended to clarify some of the problems presented 
by steady-flow systems, and to show how their 
behaviour can be investigated and quantitatively 
specified. Some of the concepts and mathematical 
expressions which appear below have already been 
used by GiLLILanD and Mason [8], [9] in the course 
of a study of flow through fluidised beds. 


F-DIAGRAMS AND AGE-DISTRIBUTION FUNCTIONS 


The volume of the vessel occupied by the fluid is 
V, and the volumetric rate of inflow and outflow of 
fluid is assumed constant and equal to v. Suppose 
some property of the inflowing fluid undergoes a 
sudden change from one steady value to another ; 
for instance, let the colour change from white to 
red. Call the fraction of red material in the outflow 
at time @ later be F (6). The plot of F (8) vs. v 0/V 


Fig. 1. F-diagrams. 
(a) Piston flow; (b) Piston flow with (c) Complete mixing; (d) Dead water. 


some longitudinal mixing; 


will be called an “ F-diagram.”’ Figs. 1 (a)-(d) show 
F-diagrams for some representative types of system. 

Perfect piston flow (1(a)) will never occur with 
Newtonian fluids; there will always be some longi- 
tudinal mixing, due to viscous effects and molecular 
or eddy-diffusion. 1 (b) illustrates the departure from 
piston flow caused by restricted longitudinal mixing. 
1 (c) is the diagram for perfect mixing ; the equation 
of the curve is easily shown to be : 


F (6) =1 — (1) 


It cuts the ordinate v#/V =1 at (1 — *), and its 


initial slope is unity. 1(d) shows a diagram of the 
type to be expected when there is a good deal of 
““ dead water ” in the system ; a considerable fraction 
of the fluid is trapped in eddies, and spends much 
more than the average length of time in the vessel, 
while most of the flow takes place through a restricted 
channel. 

The F-diagram of a system is easy to obtain (for 
instance, by injecting tracer materials into the entering 
stream), and its shape will clearly give a good deal 
of information about the behaviour of the fluid 
flowing through the vessel. As will be shown later, 
much of this information can be summarised by two 
numbers which can be derived from the F-diagram. 
The diagram also enables certain calculations to be 
made concerning the performance of the system when 
used, for instance, as a blender or reactor. 

The shape of the F-diagram depends on the relative 
times taken by various portions of the fluid to flow 
through the vessel, or in other words, on the distribu- 
tion of residence-times. Those elements of the 
material in the vessel which have been in it for a 
time @ are said to have an “age” 6, and the fraction 


of the material in the system having at any instant 
ages between @ and (@ + dé) is 1(@) dé. The fraction 
of the material having ages between @ and (@ + d@) 
at the moment of leaving the system is E(@) dé. I 
and E may be called the internal and exit age- 
distribution functions respectively. The relationship 
between E, J and F can be shown by imagining the 
ingoing stream, as before, to change from white to 
red at time 6=0. At a time @ later the balance 
sheet for the red material is : 


Entered : v6. 
Still in system: V frm dé’. 


6 
Left system : ef ae. 

=0 
(The last follows because the rate of outflow of red 
material at any time @” after the change of colour 

is Thus the conservation equation 


for red material is 


6 6 6 

v 

0 


Also, the total fraction F (@) of red material in the 
outflow at time @ is 


F (0) = de’. (3) 
oO 
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Differentiating (2) with respect to @: 
1- = (6), 
t 
and hence from equation (3) : 

1 F(6) =~ (5) 
Equation (4) gives the relationship between E (@) and 
I (6), while equation (5) shows how the latter function 
can be determined from the F-diagram. 

The same information can be obtained in a some- 
what different way, which may prove more convenient 
under some circumstances. Suppose a quantity Q 
of some substance (a radioactive tracer might be 
chosen for ease of estimation) is injected into the 
entering stream virtually instantaneously — that is, 
within a period very short compared to Vv. As 


0 
so that the area under every C-diagram is equal to 
unity. 

It has been assumed in the foregoing that the 
forms of the functions J and E do not vary with 
time. In fact, however, instability of flow may lead 
to fluctuations in the internal and exit age-distribu- 
tions, and cause observed values of F (@) to oscillate 
about some smooth curve of the kind shown in Fig. 1. 
These fluctuations will generally be random in nature, 
and by making a number of determinations of the 
F-diagram and averaging, the mean form of the 
diagram can be obtained. The problems raised by 
such fluctuations will be considered on another 
occasion ; for the moment attention will be confinea 
to cuses where they may be ignored. 


é 0 


0 


Vig. 2. C-diagrams. 


(a) Piston flow (b) 


Piston flow with 


(c) Complete mixing (d) Dead water. 


some longitudinal mixing 


before, samples are taken at various times @ after 
the injection; the concentration of the injected 
material in the exit stream at time @ is C(@). Then 
it can be shown, by methods similar to those used 
above, that 


F (0) E (6). (6) 


v dé 
VC (0) 


v6 


Tr’ which may be called 


Plots of s. 


C-diagrams, are shown in Fig. 2 for the systems 
whose F-diagrams are given in Fig. 1. Note that 


It will be noted that 


= 


v 
- Fajen (8) 


E(@)d0 =1 (9) 


@ 
(because the integrals in equations (8) and (9) repre- 
sent the total fractions of material, in the vessel 
and at the exit respectively, having ages between 


3 


| 


P. V. Danckwerts : Continuous flow systems 


0 and @o — namely, all the material in each case). 
Equation (8) implies that the area between the curve 
F (@) vs. v@/V and the line = 1 is equal to 
unity. It can be seen immediately that the two 
shaded areas are equal in Fig. 3 whatever the shape 
of the curve. 


1-A 


7 


Fig. 3. 


The average age 6, of the material leaving the 
vessel is always V/v. This can be shown as follows : 


(10) 


— | 0. 


Hence, from equation (3) 
v6, F(6),, 


It can be seen that from Fig. 8 that the last integral 
is equal to 1. 

The average age of the material in the system at 
any time is @,, where 


Integrating by parts we have 


v 

| @ d F 
F=0 

The integral in (12) or (18) can be evaluated graphic- 

ally with the help of the F-diagram. 


é, = (18) 


Hoip-Back SEGREGATION 

The systems illustrated in Figs. 1(a)-(d) show in 
progressively greater degree a quality which will be 
called “ hold-back.” In a system which displays 
hold-back, some elements of fluid spend more, others 
less, than the average time, V’/v, in the vessel. In a 
system with piston-flow (1(a)), there is no hold- 
back. Quantitatively, the magnitude of the hold-back 
will be denoted by H, and can be conveniently defined 
as the area under the F-diagram between v@/V = 0 
and v@/V = 1 (area A in Fig. 38) : 


Vie 
v 
H = | 


(14) 


H varies from 0 for piston-flow to values approaching 
1 when most of the space in the vessel is dead water. 
For the completely mixed vessel (Fig. 1 (c)), H = 1/e. 
The magnitude of H is a measure of the deviation from 
piston-flow. Its significance is most easily visualised 
as follows: if the colour of the inflowing stream 
changes suddenly from white to red, H is equal to 
the fraction of the vessel which will still be occupied 
by white fluid after a volume of red fluid equal to 
the volume of the vessel has flowed in. 

The hold-back is of importance in connection with 
chemical reactors and the successive flow of fluids 
through pipes and other vessels. In the former case 
it will usually be true to say that if the reactants are 
adequately mixed before or shortly after entering the 
vessel, the output of the reactor will be greatest when 
the hold-back is least (although the reverse may be 
true for autocatalytic reactions). A high hold-back 
means that much of the volume of the reactor is 
occupied by material which has already undergone 
reaction, while much reactant passes rapidly through 
the vessel by a “ short-circuit "’ route. The magnitude 
of the hold-back gives only a general idea of the 
behaviour of the system, however; the whole F- 
diagram (and possibly other information) is needed 
if the actual performance of the reactor is to be 
calculated. This subject is discussed in a later 
section. 

An indication of the efficiency of mixing in a system 
can be given by a single quantity, S, which will be 
called the “ segregation,” and can be derived from 
the F-diagram as follows. If one superimposes the 
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F-diagrams for a perfectly-mixed system (equation 
(1)) and for an imperfectly mixed system such as that 
represented in Fig. 1(b), the result will be similar 
to Fig. 4. 


Fig. 4. 


The degree of departure of the system from perfect 
mixing will be indicated by the size of the shaded 
area (4, + A,). However, since the area between 


each curve, and the line F (@) = 1 is the same, the 
two shaded regions have the same area (4, = Ay). 
It is thus convenient to define the segregation S 
quantitatively as the area A, between the F-diagram 
of the system and the curve F (@) = 1 — e/”, up 
to the point (at @ = T) where the curves cross; it 
is thus equal to half the total shaded area. 

When there is dead water in the system the positions 
of the curves may be inverted, as in Fig. 5. 


The segregation is then given the negative value 
— A, to indicate the nature of the departure from 
perfect mixing. In a case such as that illustrated in 
Fig. 1 (b), the two curves may cut twicc, as indicated 
in Fig. 6. 


Fig. 6. 


In such cases .1, + A, = Ay, and the segregation 
is equal to — A,. S varies from + 1/e for piston- 
flow to values approaching — 1 when most of the 
space in the system is dead water. 

The sense in which the word “ mixing ” has been 
used here should be clearly understood. We are not 
concerned with homogeneity, but with the age- 
distribution of material in the vessel and in the out- 
going stream. If these distributions are the same, 
the system is said to be perfectly mixed, although 
there might, in the experiment described, be a 
perceptible non-uniformity of colour. On the other 
hand, when a fluid flows through a pipe with a high 
degree of turbulence, the segregation as defined above 
may approach that for piston-flow, although mixing 
in directions at right-angles to the axis of flow may 
be very efficient, ensuring that at any given moment 
the composition will be uniform across each cross- 
section and in the outflowing fluid. 

As pointed out by Gitiianp and Mason [9], 
although the F-diagram for a fluidised bed may be 
quite close to that for perfect mixing, other tests 
show that longitudinal mixing is in fact not very 
efficient, and moreover a considerable proportion of 
the gas passes through the bed in the form of bubbles. 

Under some circumstances, when the fluid concerned 
is a gas, molecular diffusion between neighbouring 
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elements of fluid may play a material part in deter- 
mining the shape of the F-diagram. In liquids, how- 
ever, diffusion coefficients are so small that only 
hydrodynamic effects will usually be of significance. 


turovcn Beps or Soiips 

A very large number of chemical engineering opera- 
tions involve the ‘flow of a fluid through beds of 
stationary solid particles. Calculations are usually 
based on the assumption that piston-type flow occurs 
in such systems. This is certainly not exactly true, 
since longitudinal mixing must take place. Many 
cases, however, are probably amenable to the mathe- 
matical treatment which follows. 


Fig. 7. 


Consider a packed tubular vessel of length L, 
through which fluid flows with a mean axial velocity 
u. At time @ = 0 the colour of the fluid changes 
from white to red. Now if the flow were of piston-type, 
the plane boundary between red and white would 
move down the tube with velocity u. We shall 
denote this imaginary plane by z = 0 and use it as 
the origin of a frame of reference moving down the 
tube with uniform velocity u, so that at time @ the 
plane z = 0 is distant u@ from the entry, and the 
z-coordinates of the ends of the tube are (— u@) and 
(L — u@) respectively. (See Fig. 7). 

The longitudinal velocity of any element of fluid 
relative to the plane z = 0 will fluctuate irregularly. 
At times the element will be close to a solid surface, 
and viscous forces will slow it down, while at other 
times it will be near the centre of a channel and 
moving at a velocity greater than the mean; part 
of the time may be spent in regions of turbulence 
where the element will undergo rapid and irregular 
fluctuations in velocity. If the packing is quite 
randomly arranged, without any channelling, each 


element of fluid will travel at the same average 
velocity, and will experience fluctuations of the same 
average magnitude and frequency. It is clear that the 
ordinary “ random walk " theory [10] can be applied 
to such behaviour; the result will be redistribution 
of red and white material according to the laws of 
diffusion. That is, if ¢ is the mean concentration 
(volume fraction) of red material at a plane 7 at 
time @, 


(15) 


D is a “ diffusivity’ which must be determined 
empirically ; it will presumably depend on the 
viscosity, density and velocity of the liquid, and on 
the size and shape of the packing. The following 
boundary conditions will represent the facts with 
sufficient precision under many circumstances : 


c@=0,7>0, 

e=1,2<0, @=0 

c=0,2=0, 

c=lr=—o, @>0 
(The above conditions hold if the concentration of 
white material at the entry falls virtually to zero at a 
time very much less than L/u. It can be shown that 
this will occur if 4D/Lu <1-that is, provided 
mixing is not too effective or the tube too short). 
Under these conditions the solution to (15) is 


values may be found in tables and are shown in 
Fig. 8). Hence 


(16) 


(17) 


(18) 


L—uwu 
2F (0) =1 —erf 2\/ Dé 


and since v/V = u/L (V being the void volume of 
the vessel) : 


1 — vb, 


The F-diagram will thus be determined entirely by 


(19) 
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the value of D/Lu. It will have the general shape 
shown in Figs. 1 (b) and 9. 
If we let 


erf (z) =1 —2F (20) 


equation (19) becomes, on rearrangement, 
1—v0/V D 
Lu 
(Values of = corresponding to given values of F can 
be found with the help of Fig. 8). 


(21) 


6 


2 


Fig. 8. Erf (a) = 


For a given flow-rate and depth of packing, the 
left-hand side of equation (21) should be constant if 
the diffusion equation (15) is obeyed. However 
equation (21) is not very suitable for testing the 
diffusion hypothesis or for determining the value of 
D since the expression on the left-hand side is very 
sensitive to small errors in F over the range of greatest 
importance. The value of D is probably more easily 
determined from the slope of the F-diagram at 
vé/V = 1. Differentiating (19) with respect to 
and putting v@/V = 1, we find 

dF ral Ja 

D/Lu can be calculated from this expression, the 
F-diagram calculated from equation (19), and the 
results compared with the experimental curve. D 
should, of course, be independent of L for a given u. 


(22) 


The hold-back in the system is equal to »/DV/zv, 
or DL 

Fig. 9 shows the F-diagram obtained when water 
flowed through a bed of j-in. Raschig rings (tube 
diameter 4.8 cm; L = 140 cm ; fractional free volume 
0-62; volumetric flow-rate 4-5 cc/sec). The curve 
calculated from equation (19) with D/Lu = 0-018 is 
compared with the experimental points and shown 
to agree well. 


| 


08 70 
Fig. 9. F-—diagram calculated for flow through packed 
column compared with experimental points. 
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The liquid running over the packing in an absorption 
tower is presumably subject to the same kind of 
mixing, although channelling may prevent equation 
(15) from being very closely obeyed. The quantity V 
in this case is equal to the total liquid hold-up. 

The effect of longitudinal mixing on the performance 
of tubular reactors is discussed later. 


FLow 1s Pipes 

Hold-back will always arise when a true fluid flows 
through a pipe. In the first place, there is a variation 
in velocity from the axis to the wall of the pipe, so 
that the central “core” of fluid moves with a velocity 
greater than the mean, while the fluid near the wall 
lags behind. This effect will be most marked when 
flow is laminar; in turbulent flow the velocity is 
more uniform across the pipe, and the hold-back is 
less, although eddy-diffusion contributes to the 
longitudinal mixing. In some circumstances molecular 
diffusion may also contribute appreciably to the 
hold-back. It will be convenient to discuss laminar 
and turbulent flow separately. 
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When viscous liquids such as heavy petroleum oil 
or rayon dope are pumped through pipes the flow 
may be laminar. If the flow is changed, say, from 
one type of petroleum to another, or from normal to 
pigmented dope, the mean composition of the liquid 
leaving the pipe will change only gradually from one 
pure component to the other, and a large volume of 
mixed liquid will be discharged. If certain simplifying 
assumptions are made, the form of the F-diagram 
is easily calculated. (Similar calculations have been 
made by Boswortn [2], and by Fow.er and Brown 
(7).) 

In the first place, entrance-effects will be ignored 
(as in the derivation of Poiseuille’s formula), and the 
liquid will be assumed to be everywhere in unacceler- 
ated laminar flow with a parabolic velocity-distri- 
bution; provided the ratio of the length to the 
diameter of the pipe is sufficiently large (> 0-06 Re) 
entrance effects will in fact be unimportant. Secondly, 
it will be assumed that the two liquids have the 
same density and viscosity, and behave as Newtonian 
fluids. Thirdly, molecular diffusion will be ignored. 

If & is the mean velocity and u the velocity at a 
distance r from the axis of a pipe of radius R, we 
have the well-known relationship 


w= 5), (28) 


Ly 
L /, 


Fig. 10. 


The follow-up liquid will thus first appear at the 
exit of a pipe of length L after a time L/2u; after a 
greater time, @, the situation will be as shown in 
Fig. 10, with a core of follow-up liquid of diameter 
2R 1/1 — L/2U8 at the exit. Allowing for the varia- 
tion of velocity from the centre to the edge of the 
core, the fraction of follow-up liquid in the stream 
leaving the pipe is seen to be 

RV 


F (6) = “ Qnr.dr, (24) 


= 1 — 6, > 4 
=1 — 00/V 


(25) 


os 
Wy 


Fig. 11. F—diagram for viscous flow in pipe. 


The resulting F-diagram is shown in Fig. 11. Its shape 
is independent of the length and diameter of the pipe, 
and of the viscosity and velocity of the fluid, (provided 
of course that flow is laminar). The hold-back is }. 
Boswortn [2] has calculated the effect of molecular 
diffusion on systems of this kind (his F, is equivalent 
to the present author’s E (@)). He shows that it can 
be ignored providing the following conditions are 


fulfilled : 
R > 18)/D,, Lju 
L>65 x 


(where D,, is the molecular diffusivity), and points 
out that in the case of gases turbulence usually occurs 
while R is too small for the first condition to be 
satisfied. With liquids, however, which have much 
higher Schmidt numbers (u/pD,), there may be a 
considerable range of pipe-diameters for which flow 
is laminar and the first condition is also satisfied. 
Boswortu [3] has also derived expressions for E (@) 
(F, in his nomenclature) for turbulent flow in pipes, 
starting from simplified expressions for the values of 
the axial velocity and eddy-diffusivity at various 
distances from the pipe wall. The following expression 
for F (6) can be derived from Boswortn’s formula : 


= 1 — 2a” x 


— B(2n + 1) (a 
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x V 
n 1, 5B 
x — pn — p/2) + (—5+ x 
2n Bn 4n 
B =(n — 1)? R/0.32 nL, 
n = 7, 2,000 — Re <— 100,000, 
n = 8, Re > 100,000, 


R being the radius and 1 the length of the pipe. 


10 
Fig. 12. Experimental (A) and calculated (B) F—diagrams. 
Water in circular tube, Re = 6940, R/L = 63 = 10+. 


1S 


Comparison can be made with the experimental 
results of FowLer and Brown [7]. Fig. 12 shows an 
experimental F-diagram compared with that cal- 
culated from Boswoxrtn’s formula. In Fig. 13 a 
comparison is made between calculated and observed 


values of — for various Reynolds 


numbers in a given pipe (i.e. F has the value 0-8 
at 6,,, and the value 0-2 at 6,,). It will be seen 
that Boswortn’s formula predicts (at least under the 
conditions chosen for the comparison) a much greater 


degree of longitudinal mixing than is actually ob- 
served, and also that it fails to predict the marked 
effect of velocity of flow on the F-diagram for a 


given pipe. 


2000 4000 8009 16090 32000 
Re 
Fig. 13. Experimental (A) and calculated (B) values of 
v/V — % 9) RIL = 63 x 16+. 


Reactors 


The performance of a steady-flow reactor can be 
calculated from the F-diagram provided the reaction 
is first-order, or pseudo-first-order, with a velocity- 
constant which does not vary from place to place 
in the reactor. For instance, if an effluent containing 
a radioactive element is allowed to flow through a 
reservoir, the concentration of the element at the 
exit can be calculated if its radioactive decay-constant 
k and the F-diagram of the system are known. 

Consider those elements of fluid which have ages 
between @ and (@ + d@) at the exit. A fraction 
(1 — e~**) of the reactant originally contained in 
these elements will have reacted during its passage 
through the vessel. Hence the total fraction (1 —f) 
of reactant which reacts during its passage through 
the vessel is 


(27) 


1 -s=fa E (6) 


From equation (3) 
E (6) d0 =d F (6). 
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1 fa dF (6) 


=} — few ar (a 
6=0 


f can be evaluated by graphical integration if the 
F-diagram is available. 

In the special case, previously discussed, of a 
tubular packed vessel, the expression for F given in 
equation (19) is unsuitable for use in equation (29) 
because of the approximations which have been made. 
Instead we set up the differential equation for a 
tubular reactor with longitudinal diffusion as well as 
flow (changes in volume are assumed not to occur, so 
that the mean longitudinal velocity, u, is the same 
at all cross-sections). The equation is easily shown, 
by making a balance on a section of differential length 
dy, to be 


(30) 


wherecis the concentration of reactant at cross-section 
y. The concentration of reactant in the entering 
stream is c*; owing to diffusion the concentration just 
within the entrance of the reactor, at y = 0, is less 
than c*. The following boundary condition expresses 
the fact that the rate at which reactant is fed to the 
reactor is equal to the rate at which it crosses plane 
y = 0 by combined flow and diffusion. 


uc* = uc -- D.de/dy,y =0 (31) 
At the outlet (y = L) 
ufe* = ue — D. de/dy,y = L, 


fc* being the concentration of the exit stream. Now 
if de/dy were negative, the concentration in the exit 
stream would be greater than that at the end of the 
packing. If de/dy were positive the concentration 
would pass through a minimum somewhere in the 
reactor and then rise towards the downstream end. 
Intuition suggests that neither of these situations 
can arise, so that the boundary condition must be : 


de 
(32) 


The solution to equation (30) with these two boundary 


2(1 +a)exp (L —y)—2(1 —a)exp ly —L) 
x 


(1 +a)? exp — (1—a)*® exp D 


where 
a=7/1 + 


The value of (1 — c/c*) at the exit (y = L) is equal 
to (1 —/f), the fraction reacted. Hence 


1-f=1- 
; 
(1 +a}exp (1—a)—(1 —a)exp “(1 +a) 


As D +0 this expression tends to 


1—f=1—exp = (35) 


which is the well-known solution for “ piston-flow.” 
As D - ©, we find 
kL 
36 
1-f=— (36) 
which is the solution for complete mixing. For very 
small values of D 


1 —(1 DL (37) 


From this it can be seen that the effect of diffusion 
is to decrease the fractional conversion compared to 
that for piston-flow, but that the effect will be negli- 
gible provided k? DL/u® <1, or [Inf]? D/Lu <1. 
Thus for specified operating conditions it is possible 
to decide, if an approximate value of D is available, 
whether neglect of longitudinal diffusion will lead to 
serious error in calculating the output of the reactor. 

The foregoing treatment of reactors will seldom be 
applicable to practical problems except as a rough 
guide. In the first place, the heat of reaction often 
gives rise to temperature-gradients, and hence point- 
to-point variations in the value of the reaction- 
velocity constant, so that the chance of a molecule 
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reacting depends on its path through the reactor, as 
well as its residence-time. Secondly, if the reaction 
is of order other than first, the chance of a given 
molecule reacting depends on the molecules which 
it encounters in its passage through the reactor ; the 
nature of these encounters is largely determined by 
diffusional processes caused by point-to-point varia- 
tions in composition in the fluid, which cannot be 
deduced from the F-diagram. Second-order reactions 
between imperfectly mixed fluids have been consi- 
dered by the author [5]. Densien [6] has considered 
second-order reactions taking place in a fluid in 
laminar flow in a pipe, under conditions such that 
molecular diffusion may be ignored. 


BLENDERS 


A continuous-flow blender is a mixing vessel into 
which flows a stream of material of continuously- 
varying composition. In the vessel elements of 
material which have entered at different times are 
mixed, so that the outilowing stream shows less 
variation in composition than the input. The case 
of a perfectly-mixed blender has been considered in 
detail elsewhere [4], with special reference to the flow 
of town gas through a gas-holder, and consequent 
smoothing-out of fluctuations in the calorific value. 
Beaupry [1] has considered the special case of a 
perfectly-mixed blender with a feed consisting of 
batches of material of finite volume. A_ blender 
which displays segregation will be less effective in 
reducing fluctuations than one in which perfect 
mixing occurs. 

Suppose the concentration of some component in 
the entering stream has a concentration ¢; which 
fluctuates with time. Its value at any time is 
(@ + 8,), where @ is the mean value of ¢; which is 
assumed to show no trend. A convenient measure 
of magnitude of the fluctuations in ¢; is its standard 
deviation : 

o? = 82. (38) 
a; is assumed to show no trend. The concentration 
Cy = (@ 4- 8) of the outgoing stream has a standard 
deviation o,. The problem is to find the ratio a, 9;. 
The rates of inflow and outflow are assumed constant 
and equal to v. The concentration c, (t) in the outflow 
at time ¢ is then given by 


(t) = {t — 0} E(0) dé, (39) 
6=0 


where c;{t — 6} is the concentration in the inflow 
at time (t — 6). Hence 
@ 
5, (t) = fa {t — 0} dé 


(40) 


(3, (0)? = [fact ~ 6} ae). (41) 
6=0 


However, it is generally true that 
@o o 
fw) ay| | | Swfly +r) dydr, (42) 
0 y=0r=0 
where f(y) is any function of y. Applying this to 
equation (42) : 
o @ 
(t)]? = 2 {t — 0} x 


x 8 {t — —r} E(0) +r) dédr. 
Averaging this with respect to t 


og = 83 =2 — 0 x 


<x E(@) E(@ +r) 


(44) 


Now the quantity 
R(r) = 


(the numerator being averaged with respect to t’) is 
known as the autocorrelation coefficient or serial 
correlation coefficient of c; for a time-interval r. Sub- 
stituting (45) in (44) we have 


(45) 


R (r) E(0) +7) d@dr. (46) 


6=0r=0 


— 
o,2/07 = 2 


R (r) may be found as a function of r from a repre- 
sentative record of c; vs. t in various ways which have 
been discussed elsewhere (4). (The relationship 
between R (r) and r indicates whether the fluctuations 
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in ¢; are on the whole rapid or slow, and whether 
they display any regular periodicity). Assuming that 
R (r) is known for a sufficient range of values of r, in 
order to evaluate the integral in (46) it is necessary 
to know E(@) as a function of 6. The necessary 
information is contained in the F-diagram of the 
system, and the following is perhaps the simplest 
way of doing the calculation. First fit some simple 
function of @ (such as a power-series) to the F-diagram 
up to a value 6,, at which F (6,) is very close to 1. 
Differentiate this function with respect to @. From 
equation (3) we see that 

dF (@) 

E (8), (47) 
so that E(@) is now a known function of @. For a 
given value of 6, E(@ +r) is the value obtained by 
replacing @ by (@ +r). It will now be possible to 
evaluate the integral 


a 
I(r) = r) dé (48) 
6=0 
algebraically for a number of values of r. The integral 
o,?/02 = 2 fe (r) 1 (r) dr (49) 


r=0 


can now be evaluated graphically. (In practice, of 
course, the upper limit will be taken at some value 
of r for which the integrand has become vanishingly 
small). 


Tue Use or Mopets ror Prepicrinc F-pIAGRAMS 


Under certain circumstances the F-diagram of a 
large system can be predicted with confidence from 
that for a model. The conditions to be fulfilled are 
(a) The model must be geometrically similar to the 
system ; (b) The Reynolds number must be the same 
in the model and the system; (c) Gravity-waves, 
density-differences, surface tension, and other in- 
fluences apart from inertia and viscosity, must be 
unimportant in determining the behaviour of the 
fluid in both model and system. 

In the simple case of a fluid flowing through a tank 
or a fixed bed of solids, the same fluid being used in 


the model as in the system, equality of the Reynolds 
numbers requires that 


+= 50 
(50) 
where v,, v, are the volumetric flow-rates and L,, L, 
the linear dimensions of the model and system 
respectively. That is, in a 1/10 scale model the 
volumetric flow-rate must be 1/10 that in the full- 
scale system. If the system incorporates a rotating 
stirrer, the angular velocities, w, of the stirrers in 
model and system are related by 


In the 1/10 scale model the stirrer must rotate 100 
times as fast as in the system. 

If these conditions are fulfilled the flow-patterns 
of the fluid will be geometrically similar in model 
and system, hence both will have the same F-diagram. 

The simple type of model experiment described 
here cannot be applied to systems in which the fluid 
is of non-uniform density or viscosity owing to varia- 
tions in temperature or composition from point to 
point, or to fluidised beds. 
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NOTE ADDED IN Proor 
Flow in Pipelines—Huwrt and Kent [11] have recently 
published an account of experiments on petroleum flowing in 
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a 182-mile long 10-inch diameter pipe line. A_ radioactive 
tracer was used to obtain (in effect) C-diagrams for various 
lengths of pipe. With a Reynolds Number of 20,000-30,000 
the diagrams obtained were close to those which would be 
expected if the * diffusion equation (15) were followed, D 
having « value of 1-3 sq. ft. /sec. Further experimental work 
on mixing in pipe-lines is referred to by the authors [12], [13]. 

Flow in Packed Beds—Lariwvus and Amunpson [14], 
have analysed the consequences of longitudinal mixing in 
adsorption columns. Breannano and Wituetm [15] have 
measured the apparent radial diffusivity in fluids flowing 
through beds of solids, It seems likely that the longitudinal 
diffusivity D (equation [15]) will have a different value. 
Anruur et al. [16] and Moraces et al. [17] have shown that 
the fluid velocity differs markedly from the mean in the 
neighbourhood of the wall of the column. 


NOTATION, 


A = area under curve 
¢ = concentration 
C(@) = concentration at exit of vessel at time @. 
c* = concentration of stream entering reactor 
D = apparent diffusivity 
D,, = molecular diffusivity 
E (@) = distribution-function for residence-times 6 


ert) = 


F (8) = fraction of material in outflow which has been in 
system for a time less than 6 
J = fraction of reactant unreacted 
H = hold-back (equation (14)) 
distribution-function for “ages,” 
system 
= first-order reaction-velocity constant 
length of tube or reactor 
parameter in equation (26) 
= quantity of tracer material iniected 
radius of pipe 
serial correlation-coeflicient of c; for time-interval r 
radial distance from tube-axis ; or time-interval 
Reynolds number 
segregation (Figs. 4, 5, 6) 
time 
velocity of flow 
mean velocity of flow in open pipe 
volume of vessel 
= volumetric flow-rate 
= distance referred to coordinates moving with velo- 
city u 
distance referred to stationary coordinates 
= defined by equation (20) 
parameter in equation (26) 
= parameter in equation (26) 
= instantaneous deviations of ¢;, cg from ¢ 
@ = time-interval, residence-time or “age” 
= standard deviations of ¢;, Co from 
w = angular velocity 


6, of material in 
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30°C, 40°C, 50°C and 60°C, 


Extraction of acetic acid from water. 1—benzene-acetic acid-water 


F. H. Garner, 8S. R. M. Exuis and U. N. G. Roy 
The Chemical Engineering Dept., The University, Birmingham, 15 
( Received 4 November, 1952) 


Summary — Liquid-liquid equilibrium data are given for the system benzene-acetic acid-water at 
The experimental distribution coeflicients are in close agreement with those 


predicted from binary vapour-liquid equilibrium data; Other papers are to follow giving liquid-liquid 


and vapour-liquid equilibrium data for the system ethyl acetate-acetic acid-water. 
Résumé— Présentation de données expérimentales sur l'équilibre liquide-liquide dans le systéme benzéne- 


acide acétique—eau a 30°C, 40°C, 50°C et 60°C, Les coeflicients de partage expérimentaux sont en accord 
atroit avec lex valeurs calculées & partir des données d'équilibre liquide-vapeur des systémes binaires. Les 
résultats d'équilibre liquide-liquide ct liquide-vapeur pour le systéme acétate d'éthyle-acide acétique-eau, 


seront publiés ultérieurement. 


Introduction 


The strong affinity of acetic acid for water from which 
it cannot be easily separated has led to a considerable 
amount of research into the problem of devising 
separation methods that can be economically used by 
industry. Large quantities of acetic acid in aqueous 
solution are produced in the form of pyroligneous 
acid from the wood distillation industry, and, in still 
weaker form, from fermentation processes. Large 
quantities of glacial acetic acid used in the cellulose 
acetate industry undergo dilution and require con- 
centration before they can be reused. 

Goerinc [1] in 1883 applied liquid-liquid extraction 
to the concentration of acetic acid from aqueous 
solutions. The solvent used was ethyl acetate and 
Gorrinc drew attention to the fact that change in 
temperature had little effect on the extraction 
process. Ricarp and Guixor [2] and Dreyrus and 
Haxry make use of a second liquid solvent such as 
a hydrocarbon to reduce the solubility for water of 
the solvent layer. Ornmer, Wuire and Truecer 
provide a large amount of information on the use of 
different solvents for extracting acetic acid from 
water. 

B.1.0.S. Report 1051 mentions the use of a solvent 
85% ethyl acetate, 15% benzene. An interesting 
observation is the power of traces of methylene 
chloride in aiding the separation of ethyl acetate 
from water [4]. 

Ternary liquid-liquid equilibrium data has been 
published for systems involving acetic acid, water, 
and a large number of solvents including chloroform, 
isopropyl! ether, methyl isobuty] ketone, benzene, ethyl 
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acetate, etc. Treysa. [9] has shown a large discre- 
pancy between the predicted and experimental 
results for the system benzene-acetic acid-water. 

The purpose of this first paper is to present new 
experimental data for the system benzene-acetic 
acid-water at a range of temperatures and to compare 
the experimental and predicted distribution coefli- 
cients. 


Determination of Solubility Curve 

The apparatus used to determine the solubility curve 
was essentially a thin walled glass tube, fitted with a 
cork through which is introduced a glass stirrer and 
thermometer graduated in tenths of a degree centi- 
grade. Through another hole in the cork a burette 
is introduced so that the tip of the burette is just 
projecting beyond the cork. The apparatus was 
placed in a thermostat controlled to + 0-1°C with 
illuminated windows to view the experiment. Weighed 
miscible mixtures were placed in the test tube and 
titrated to turbidity with a third component. 


Determination of Tie Lines 


For acetic acid-benzene-water at 30°C and 40°C the 
apparatus described by Hunrex and Nasu [5] was 
inserted in a large thermostatic bath and used for 
the determination of tie lines. At temperatures of 
50°C and 60°C the liquid mixture was weighed in a 
stoppered flask and introduced into the constant 
temperature bath. The flask was agitated thoroughly 
every 15 minutes for 3-4 hours. After standing over 
night samples were then drawn off from the two 
layers. 
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Purity of Reagents 

Water—Distilled water was treated with solid barium 
hydroxide to remove any carbon dioxide. 

Acetic acid—Analar acetic acid was fractionated with 
ethylene dichloride and then stored in the absence 
of moisture. B.P. 118-15°C. 

Benzene—The purification of benzene was effected in 
a thirty plate column at a reflux ratio of 25:1. The 
boiling point at 760 mm was 80-1°C and the refractive 
index D® 1.5011. 


Analysis of Acetic Acid 

A volumetric method of analysis was used where a 
dilute solution of acetic acid was titrated with standard 
carbonate free barium hydroxide. Phenolphthalein 
was used as the indicator. 


Results 

Solubility curve and tie line data are given on Tables 
1 and 2. Fig. 1 shows the solubility and tie line data 
at 30°C represented on a triangular graph. 


Acetic acid 


Woler 
Fig. 1 


Discussion of Results 


A comparison of the results of Tables 1 and 2 shows 
that the spread of the solubility curve decreases only 
slightly with temperature. The distribution coeflicient 
shows only a small variation with temperature when 
the concentration of acetic acid in the water layer 
is less than 50-0% by weight. Above this concentra- 
tion of acetic acid in water the distribution coeflicient 
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The tie line data is satisfactorily correlated by 
Bacuman [6] and Orumer-Tostas [7] plots. 

On Fig. 2 a plot has been made of the mol fraction 
of acetic acid in the benzene layer against the mol 
fraction of acetic acid in the water layer. The resulting 
distribution coefficients are slightly lower than those 
reported by Hanp [8]. Treysat [9] states that 
predicted distribution coefficients from the available 
vapour-liquid equilibrium data were considerably 


Table 1. Acetic Acid-Water-Benzene—Solubility Data 
Acetic | Acetic 
°c | acia | Water) Benzene | acid | Water | Benzene 
| gm gm gm 
30-0, 29-3 | 69-4 0-4 40-7 | 57-3 2-0 
39-6 | 59-2 12 498 | 47-5 2-7 
48-7 | 49-1 2-2 57-7 | 36-2 61 
57-2 | 39-2 3-6 64-0 | 25-2 10-8 
642 | 29-0 6-8 592 | 139 26-9 
66-7 | 183 | 150 47-7 80 44-3 
620 | 126 | 35-4 35-6 3-6 60-8 
482 | 63] 45-5 41-0 48 54-2 
426 | 47 | 52-7 48-9 7-9 43-2 
46-7 | 58] 47-5 555 | 11-3 33-2 
520 | 80) 400 60-9 | 15-0 241 
543 | 90 36-7 63-5 | 21-1 15-4 
578 | 104! 31-8 63-0 | 25-9 11-1 
61-5 | 11-3 | 272 60-1 | 326 7-3 
637 | 139) 22-4 25-3 3-3 71-4 
65-7 | 163 | 180 260 | 71-7 2-3 
66-0 | 245 95 
62-7 | 31-7 5-6 
50-2 | 48-0 18 
40) 849 
202 | 25 673 
400/425 | 49!) 526 60-0]] 406 5-6 53-8 
| 414 48-2 9-1 42-7 
S74) 110) 316 544 | 132 32-4 
628 | 148 | 22-4 599 | 16-3 23-8 
648 | 214) 138 622 | 231 14-7 
642 | 26-6 9-3 606 | 28-8 10-6 
50-8 | 348 6-4 546 | 39-2 6-2 
40-8 | 57-5 1-7 40-1 | 56-5 3-4 
500 | 47-5 2-5 495 46-0 45 
58-5 | 36-6 49 57-1 | 35-6 73 
| 647 | 260 0-3 61D | 25-4 13-0 
(615 | 145! 260 535-8 | 13-3 30-9 
39-4 | 38 | 568 44-0 73 48-5 
519 | 88 | 393 30-4 3-6 66-0 
285 | 25!) 69-0 22-8 3-3 73-9 
28-1 | 703 1-6 23-1 | 2-8 
152 | 1-7 | 73-1 


varies more widely with temperature. 
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higher than the experimental values of Hanp [8]. It 
was thought these discrepancies might be due to 


inaccurate vapour-liquid equilibrium data for the 
benzene-acetic acid system and thus in this paper 


new equilibrium data has been used to predict the 
distribution coefficients. 


o6 ¥ 

Treybol [9] 

Experimento/ 
Predicted 


Prediction from Binary Vapour-Liquid Equilibria 
At equilibrium the activities of acetic acid in the 
two phases are equal 


= 72 
where z, = mole fraction of acetic acid in benzene. 
1, = mole fraction of acetic acid in water. 
¥; Yg = activity coefficients of acetic acid in 


Mol. Fraction acetic acid in benzene 


ak benzene and water respectively. 
If the partition or distribution coefficient is defined 
VOL, 
0 01 Of O03 O08 OS O08 as K then 2 
rachon acer 
Mol acetic acid in water K % 1953 
Fig. 2 
Table 2. Acetic Acid-Water-Benzene Tie Line Data 
Overall composition | % acid in % acid in Plait point 
water rich benzene rish 
% acid % benzene % water phase phase % acid % benzene % water 
30-0 18-6 28 | 23-4 23 46-8 46.8 6-4 
36-2 29-0 48 488 8-7 
578 17-7 245 646 20-1 
370 30-0 13-0 73-5 28-0 
40-0 113 57-0 31-7 23-5 24 471 45-5 
26-6 41-7 81-7 42:5 67 
35-0 30-5 345 47-5 
0-0 34:5 25°5 56-0 12-4 
48:7 41-0 10-3 640 22.0 
33.0 28.6 18-4 64-1 22:5 
50-0 6-25 62 47-5 0-62 44-2 8-0 
16-4 39-9 43-7 25-6 29 
36-5 20-4 341 48:5 
468 247 28-5 57-2 14-6 
57-5 22-1 20-4 63-5 23-7 
60-0 58 43-4 50-8 8-8 0-56 48-2 42-9 8-7 
134 34-8 49:8 22-4 2-34 
25-3 35-8 38-9 36-5 5-7 
34-2 32-8 33-0 470 8-75 
| 510 28-9 20-1 613 22-0 


In the above table all percentages are by weight 
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Table 3 
Mole ©, in Benzene Phas Mole ©. in water phase 
-— T Experimental Predicted 
Acetic «leetic distribution distribution 
acid Benzene Water acid Benzene Water coefficient coefficient 
31-2 61-3 75 45-0 85 46-5 0-68 0-66 
235 69-6 oo 40-0 +0 560 0-60 H 0-56 
10-5 85-4 22-8 76-5 O-7 0-48 045 
95-2 17 9-0 90-8 0-38 0-38 


The activity coeflicient of the solute in ternary 
mixtures can be predicted from binary vapour-liquid 
equilibrium data by linear interpolation or by the 
use of Margules, Van Laar, or other ternary equations. 
The fact that the vapour-liquid equilibrium data is 
determined at a different temperature to the liquid- 
liquid equilibrium does not in general greatly influence 
the predicted distribution coeflicient, since the ratio 
of the activity coeflicient varies only slightly with 
temperature. 

A difliculty in this particular application of pre- 
dicting the activity coeflicient of acetic acid in the 
benzene rich layer is the possibility of association 
of the acetic acid to give the dimer. Hanp [8], 
however, has shown that the simple distribution law 
for acetic acid in benzene holds in concentrated 
solutions but breaks down in dilute solutions. He has 
therefore postulated the presence of water in the 
benzene layer breaks down the association of acetic 
acid to give the monomer. Accordingly the vapour- 
liquid equilibrium and liquid-liquid equilibrium data 
have been based on a MW of 60 for acetic acid. 

In these laboratories Peance [10] has recently 
investigated the vapour-liquid equilibrium data of the 
system benzene-acetic acid. This data is in agreement 
with that of Rosanorr and Eastey [11] when based 
on the same molecular weight for acetic acid. They 
used a molecular weight of approximately 100 for 
acetic acid, this being calculated from the vapour 
densities of acetic acid, 

We have used the vapour-liquid equilibrium data 
of Orumer and Gitmor [12] for the system acetic 
acid-water. 

Thus using the method of linear interpolation as 
proposed by and [18] it can be 


seen from Fig. 2 that the new experimental results 
are in close agreement with the predicted values. This 
agreement is perhaps best illustrated by Table 3. 

From Fig. 2 it can also be seen that there is a large 
difference between the predicted results of TreyBaL 
[9] and those given in this paper. It would appear that 
TreyBat has used vapour-liquid equilibrium data 
for the system acetic acid-benzene based on a mole- 
cular weight of approximately 100, and for the system 
acetic acid water based on a molecular weight of 
60-0. As close agreement as shown on Table 3 can 
be obtained by using a molecular weight of 100 for 
both binary systems. 

The general conclusion of this work is that the 
experimental distribution coefficient for acetic acid 
in benzene and water is slightly lower than that 
reported by Hanp [8] but is in close agreement with 
predicted distribution coeflicients calculated from 
vapour-liquid equilibrium activity coeflicients. 
Further, the distribution coefficients are relatively 
constant for temperature changes of 30°C to 60°C. 
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The absorption of chlorine from air by solution of 2-ethyl hexene-1 
in carbon tetrachloride 


G. H. Roper 
School of Chemical Engineering, N.S.W. University of Technology, Broadway, N.S.W. 


(Received 19 August 1952) 


Summary— The disc column designed by Steruens and Morris [14] is used to study the rates 
of absorption of chlorine from air into solutions of 2-ethyl-hexene-1 in carbon tetrachloride. The 
olefin bulk concentration Cp is varied from zero to 0-034 Ib, mol per ft.3, the reaction rate constant 
k, in the ratio 20 : 1 by varying the concentration of the catalyst (iodine) from zero to 0-001 Ib. mol 
per ft.4, the ratio of the bulk concentration of the olefin Cz to the concentration of dissolved 
chlorine at the interface C; from zero to 280. 

The fractional increase in the liquid film coefficient due to chemical reaction varies from 
zero to 19, and was found to be : 


ky 


(A) 


(Ce + 000005) Cp 
where kh, is the liquid film coefficient and k’, is the coefficient for physical absorption at the 
same liquor rate, and C, is the concentration of iodine, Ib, mol per ft.5. 

Additional data are presented for the gas film coefficients of the disc column, based on the 
rates of evaporation of water and of carbon tetrachloride. 
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Résumé—A [aide de la colonne a disques de Sreruens et Morais [14] lauteur étudie la vitesse 
absorption du chlore mélangé Tair par des solutions d'éthyl 2 hexéne 1 dans CC1,. La 
concentration globule Cy doléfine varie entre 0 et 0,034 Ib. mol/ft.4 (0,54 mol gr/lit.) ; la con- 
stante de réaction &, est variée dans le rapport de 20 4 1 en modifiant la concentration du 
catalyseur (iode) entre 0 et 0,001 Ib. mol ‘ft.3 (0,016 mol gr/lit.) Le rapport entre les concen- 
trations Cy et celle du chlore dissous a la couche de passage varie entre 0 et 280. 

Le coefficient déchange liquide A, du a la réaction chimique variant entre 0 et 19, on 
trouve la relation (A), of Ay est le coeflicient déchange du a labsorption physique et Cg la 


concentration en iode. 


A partir des vitesses d'évaporation pour et CCl,, Pauteur présente des données sur le 
coeflicient d'échange gazeux pour le colonne a disque. 


of chlorine in ferrous sulphate solutions and by 
Srecror and Dopce [13] for the preparation of 
carbon dioxide-free air using hydroxide solutions. 
Theoretical studies have been made of the case of 
a gaseous solute of limited solubility diffusing into a 
liquid solution. one component B of which reacts with 
the gaseous solute [4], [12]. The reaction may be 
first or second order, but must be irreversible. The 
rate of diffusion of the gaseous solute A into an ele- 


INTRODUCTION 


For physical absorption design procedures are based 
on the two-film theory of Brunner [1] and Wurrman 
[15] for mass transfer. There are no general equations 
available relating the properties of the equipment, of 
the materials used and the gas and liquor rates, so 
that each system has its own equations characterising 
its performance. When the complication of a reaction 
in the liquid phase is encountered there is no design 


method available because theoretical treatment has 
been made for special cases only, and for these cases 
the theory has not been tested. The most reliable 
method for good design of chemical absorption equip- 
ment is to determine the overall mass transfer data 
from the operation of semi-works scale towers, as has 
been done by Riccie and Terr. [10] for the absorption 


ment of the liquid film at a distance z from the 
gas-liquid interface is : 


(1) 


The rate of elimination of the gaseous solute in any 
element of thickness dr isk, C, C,A dx. The rate of 
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diffusion out of the element of the liquid film is : 


(oC, 
N,= A D4 \ + de 
dC, Cy dx (2) 


The solutions of equations (1) and (2) have been 
accomplished for the case in which the reaction is 
instantaneous. The liquid film coeflicient for chemical 
absorption is : 

where k’, is the liquid film coefficient for physical 
absorption under the same operating conditions. For 
other rates of reaction, the chemical absorption 
coefficient may be expressed as an indefinite function 
of a number of variables, viz : 


ky =f (Re,, Se,, hee, Cy» Cy Dy. D,) (4) 


Previous experimental determinations confirm the 
dependence of the rate of chemical absorption on 
some of the factors listed in equation (4). The effect 
of the Reynolds number Re, has been examined by 
a number of workers [2], [3], [10], [13], [14]. 

Jenny [8], using a stirred batch absorber, found 
that the rate of absorption of ethyl acetate in water 
was unaffected by the presence of sodium hydroxide, 
and that appreciable time was needed for the reaction 
between the ester and the alkali to go to completion. 
The reaction between methyl formate and sodium 
hydroxide in aqueous solution is fast. An increase in 
the caustic strength of the liquor in the absorber 
increases the absorption rate. The results lend 
qualitative support to the theory that an increase in 
the reaction rate constant will increase the liquid film 
coefficient. 

A number of papers report data showing that the 
concentration of the liquid phase reactant has a large 
effect upon the rate of absorption. The absorption of 
chlorine in ferrous chloride solutions was studied by 
Srepuens and Morais [14] using a disc column. The 
effect of the ferrous chloride concentration C,, cannot 
be separated from the effect of the concentration of 
dissolved chlorine at the interface, C;. They varied 
the ratio C,/C,; from zero to 1,200. The fractional 
increase in the liquid film coetlicient due to the 
reaction was found to be : 


(5) 

The rate of absorption of carbon dioxide in alkaline 
hydroxide solutions increases with hydroxide con- 
centrations at low concentrations [6], [8], [13], [16]. 
The rate of absorption of carbon dioxide in other 
alkaline solutions is inconclusive, but it appears that 
an increase in pH is accompanied by an increase in 
absorption rate [2], [8], [5], [6], [16]. 

There are no reported data showing the effect of 
variations in the diffusivity of the reactants or of 
variations in the Schmidt number upon the rate of 
chemical absorption. 


The system chlorine, 2-ethyl hexene-1, and carbon tetra- 
chloride 

The absorption of chlorine in solutions of 2-ethyl 
hexene-1 and carbon tetrachloride offers certain 
advantages. The materials used are definite com- 
pounds with known properties. By keeping the 
concentrations of the solutes low the physical pro- 
perties of the solutions are similar to those of the 
solvent. The concentrations are easily controlled and 
estimated. The specific reaction rate may be altered 
by the presence of a homogeneous catalyst (iodine) 
forming an interhalogen compound. The increase in 
the specific reaction rate is expected to be proportional 
to the concentration of the iodine intermediate 
compound, which, in turn, is expected to be pro- 
portional to the iodine concentration. 


EXPERIMENTAL ProcepuRE AND REsULTs 
Apparatus— A laboratory column with characteristics 
comparable to commercial packed towers was recently 
reported by Srerpuens and Morxts [14]. The form of 
the column used consists of 31 carbon discs, 1-45 to 
1.5 cm diameter and 0-44 to 0-46 cm thickness. The 
dises are threaded edgeways on 2 mm diameter glass 
rod and maintained at right angles by means of a 
styrene-butadiene co-polymer cement. The general 
arrangement is shown in Fig. 1. 


Gas Coerricrents. Tur Evaporation oF 


WATER AND OF CARBON TETRACHLORIDE 
In order to evaluate the liquid film coeflicients for 
the absorption of chlorine, it is necessary to correct 
the overall coeflicients K, by the relationship 


=1/Kg — Whe (6) 
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The gas film coefficients for the absorption of chlorine equation derived from measured rates of evaporation 
from air were estimated on the basis of a general | of water and carbon tetrachloride. 


Table 1. Vaporisation of water 


Water temp. temp. Air RUT, | 
Run Water Air F. F. % Vaporisation Driving | 
No. | rate velocity - rate force | kg 
Ib. (hr) (ft). ft./see. in out in | out owt | 10-4 Ub. /hr. | 17 AP 
| | | | | 
1 | | G96 72-8 60 215 | 0-151 
2 16-6 | 682 | | 691 60 259 0-35 0-224 
0370 | 685 | 684 700 69-3 64 241 0-33 0-221 
| 379 | 696 | oF 61 235 | O40 | 0-192 
5 O370) || «698 | 2-68 O27 0-299 
7 | ooo ooo v7 us 235 27 O-260 
| O370 | | v7 62 255 O-248 
12 24 Th | 72-3 39 4-06 O41 1953 
13 161 | | 72 721 728 | 95 60 288 | 0-266 
15 | 720 Ta 728 | O49 | 254 
17 S140) 0-972 736 | 89 66 0-60 0-325 
18 810 138 m4 | 721 | 70 | 739 85 $97 067 | 0-408 
19 810 217 | m2 | 736 | 39 0-66 0-468 
21 S14) 4-07 732 | Tit | 740 740 247 O69 1-08 
22 121 O51] 73606) Th 38 3-82 O53 O-218 
23 121 | 73-4 739 | 72 730 O53 | 0-421 
24 121 158 743 73-4 741 «732 89 60 10-4 0-64 
25 12) 217 73-6 | 730 743 #738 85 61 1146 0-68 O-514 
26 121 3-00 739 | 727 743 740 145 O61 O714 
27 121 $07 743 | 725 746 st 64 247 0-67 
28 | (408 0-680 736 736 746 | 743 | 315) | 0-456 
290 | (408 3-00 74:1 738 747 17-5 O44 1-20 
| (408 2:17 698 | 698 | 72-5 95 67 10-6 O31 | 103 
31 408 158 70-7 72:1 73-4 66 739) | | (0-661 
0-972 711 | 71 730 7 65 6-16 030 0-621 
3”) | Tis 720 73-3 68 121 047) «0-771 
244 | 247 716 | 716 72-9 73-6 oF 67 | O87 | 
a 2 72-1 | | 733 | 736 95 035 | O-734 
244 | 724 | 724 | 68 623 028 | 0671 
87 | 24 0-680 728 | 72:8 730 035 | 0-391 
| 3-00 73-4 7355 | %76 7 22-9 O75 0-921 
| 161 217 750 | 750 76-4 86 12-7 069 | 555 
“) 161 76-7 SS 66 74 O51 O451 
0-680 726 | 726 78-2 730 38 7-31 0-54 0-409 
300 | 705 | 00-8 77 30 13-4 O82 | 0-566 
+ 810 69:5 69:5 72 304 O77 | 154 
716 716 | 73-9 744 so 23-4 0-882 
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Table 2. 


Vaporisation of carbon tetrachloride 


CCL, temperature 
r v °F. 


Gas temperature 
F. 
mol hr. 


Ib. (hr. 


(mean) 


075 
0-070 
0-121 
0-110 

0-095 O68 
0-088 0456 


In all runs, air was passed up the column to give 
countercurrent flow. For the evaporation of water, 
air entering and leaving the column was tested with 
a wet and dry bulb thermometer from which the 
humidity was calculated. The exit gas for runs using 
carbon tetrachloride were analysed using a thermal 
conductivity bridge. The overall coefficients kg were 
calculated from the rate of evaporation, the log mean 
partial pressure difference in atmospheres, and the 
dry surface of the discs. 

Data for the evaporation of water are given in 
Table 1 and for the vaporisation of carbon tetra- 
chloride in Table 2. The data are well correlated 
by the equation : 


ke Puy = 0-15 
3600 py, 


with a standard deviation of 8%. 


[18 Seg (7) 


These data are 


1 Liquid feed vesse/ 

2 Niveau for contro/ 

J Dise column 

Adjustable lute vesse/ 
for liquid run off 

5 Air Flow meter 

6 Chlorine flowmeter 

7 Mixing vesse/ 

& Woter manometer 

9 Gas sampling points 

10 Excess liquid overflow 


Fig. 1. The general arrangement of the experimental apparatus. 


consistent with the gas film coeflicients for the 
absorption of ammonia in water reported by Sreriiens 
and Morus. They correlated their data using a gas 
velocity relative to the surface velocity of the water 
for both co-current and counter current flow. Their 
data for counter current flow have been recalculated 
and plotted in Fig. 2. It may be seen that they agree 
with the vaporisation data. This correlation could 
have been improved by using volume rate of liquid, 
but the data are insufficiently accurate to depart 
from the usual practice. 


Physical absorption of chlorine 

Air and chlorine were passed through tlowmeters and 
mixed in vessel 7 (Figs. 1). The gases were analysed 
before entering and after leaving the column. The 
liquor was analysed at the lute vessel. From the 
liquid and gas analyses the rate of absorption of 
chlorine was calculated and the run rejected if the 
calculated rate from the gas analysis differed from 
the rate calculated from the liquid analysis by more 


* Vaporisation of woter 7-9 to 408 t) 
« Absorption of NH,-I°S9 to 400 


100 200 =300 600 800 100 
Reynolds number for Cas Flow 


2000 3000 4000 


Fig. 2. The gas film cocilicients for the dise column. 


Run 
1038 OWS S46 760 2°35 
2 258 0-405 84-0 79-0 82 85 242 
3 103 10-0 85-0 73 86 
258 10-0 84-3 59-0 72 83 20-8 
103 200 84:5 71-5 78 RS O46 
258 200 84-2 740 77 
| 
IL. 
__| 
4 Gas | {+ 
out. 
Data of Stephens end Morris 
Liquid 10 | | 
| 
3 ‘5 Airin 0.08 
006 
7 | ; | 
Chlorine in 
008 
Liquid ou: 
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than 10%. The overall coefficients were calculated 
from the equation : 

Ky, N, AP,,, (3) 
From the values of A, the liquid. film coetlicients 
were calculated from equations (6) and (7). 

Previous liquid film data for this column [14] were 
correlated using the Suerwoop-HoLioway relation 
{11} in the form : 

= Re?" (9) 

Equation (9) indicates that a change of temperature 
from 65°F to 80°F should increase the liquid film 
coefficient nearly 20%. On this temperature range 
the value of the Henry's law coefficient decreases by 
20% [7]. Thus the product Hk’, should be nearly 
independent of temperature. 

The data for physical absorption of chlorine from 
air by carbon tetrachloride are given in Table 3 and 
plotted in Fig. 3. The data may be represented by : 

Hi’, (10) 
with « standard deviation of 10%. 


Chemical absorption of chlorine 

The first series of tests were made without using the 
iodine catalyst at constant liquor rate [’ —80-2 lb /(hr.) 
(ft.). The procedure was similar to that already 
described except that a small quantity of 0-1 N hydro- 
chlorice acid solution was placed in the liquid run-off 
tube. The concentration of physically dissolved 
chlorine in the liquor was calculated from the con- 
centration of chlorine in the hydrochloric acid solution 
and its distribution coefficient between the acid and 
carbon tetrachloride. This distribution coefficient for 
several temperatures was calculated from data 
reported by Perry [9] and Jnternational Criticial 
Tables [7] and was checked by making three runs 
without any olefin in the liquor. These experimental 
distribution coeflicients agreed with the calculated 
values within the accuracy desired. The total concen- 
tration of chlorine in the liquor was determined from 
the decrease in the olefin concentration after suflicient 
time was allowed for all the chlorine to react. 

There was some substitution reaction as indicated 
by the smell of hydrocholric acid gas from the liquor. 
The amount of chlorine consumed by the substitution 
is small and is neglected. The ratios of the amount of 
chlorine absorbed determined from the liquor analysis 


i 


Table 3. 
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Hk, = 0-000.” 
§ 


- Ub/(hriift) 


Fig. 8. The liquid film coefficients for the absorption of chlorine 
from air by carbon tetrachloride. 


to that determined from the gas analysis are generally 
between 0-85 and 1.05. 

With iodine in the liquor, the presence of the 
interhalogen compound prevented the use of the 
distribution method for calculating the concentration 
of physically dissolved chlorine in the liquor. The ratio 
of the physically dissolved (free) chlorine to the total 


T 


> 
= 


Tote/ absorbed chlorine 
Ss 


Physically dissolved chlorine in exit liquor 


T =60-2 \o/(hri(f) 
0-04 
\ 
002 
\ 
0-01 
0 0-7 03 


Fig. 4. The fraction of the total chlorine absorbed which is 
physically dissolved (in the absence of catalyst). 


factor (Cp, C,, C;, I’) independently ; examination of 
the results shows that Hk, increases with the olefin 


amount absorbed was esti- 

mated from Fig. 4 for liquor 
rates of 80-2 Ib./(hr.) (ft.), in 
which this ratio is plotted as a 


function of N ,/C, in the ab- 
sence of catalyst. The chlorine 
concentration ut the interface, 
C,, was obtained from the par- 


tial pressure at the interface 
required to give the observed 
rate of transfer through the 


gas film. 


Tests were made at a higher 
liquor rate of 180 Ib. /(hr.) (ft.), 


and were treated in the same 


Fractonal increase in Liquid Film Coef¥icvent due to chemical reachon 


manner. The overall coeffic- © Uncatolysed 
ient K, and the liquid film 
from equations (6). (7) and ? =802 
| | | 

Discussion oF RESULTS oat ; 1 

. 10 20 40 100 200 400 

The liquid film coefficients Hk, Concentration of olefin in Ingurd 
given in Tables 4 and 5 cannot Cencentretion of dissolved chlorine at the interface 


be correlated against any one Fig. 5. Increase in liquid film coefficient due to chemical reaction. I” = 80-2 Ib. /(hr.) (ft.). 
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a measure of the reaction rate constant 
ke in the absence of catalyst. The 
increase in « is proportional to C,. This 


result may have been anticipated because 
of the symmetry of k, and C, in the basic 
equation (2). Hence we may write : 


* 


os 
ky —1=89-2 P0005) a) 


© 
x Uncotalysed 
180 \b/(hri(ft) 


| | | 


Frachonol increase in Liquid Film Coef¥icrent due to chemical reacton 


* 


© Co =0-77 x10" mol/cuft 


k, 
C, | 


(138) 


CONCLUSIONS 


Stepnens and Morris found that the 
increase in liquid film coefficient was 
proportional to the 0.83 power of the 
ratio C,/C;. In this investigation the 


N 


10 20 40 100 


Concentration of olefin in hauid 


Concentration of dissolved chlorine at the interface 


Fig. 6. 


reaction. I” 180 Ib. /(hr.) (ft.). 


concentration, but decreases as the concentration of 
chlorine at the interface increases. It was therefore 
decided to attempt a correlation in the form used by 
Srepuens and Mornis although the latter was derived 
on the assumption of reaction rates infinitely fast 
compared with rates of diffusion, which leads theoreti- 
cally to equation (5). The fractional increase in the 
liquid-film coeflicient, by this treatment, is a function 
of the ratio between the concentration of dissolved 
reactant (C,,) and that of dissolved chlorine at the 
interface. 

The data of runs 1 to 19 (Table 4) are plotted in 
Fig. 5, and the data for runs 20 to 29 in Fig. 6. The 
data fall on parallel straight lines, having a slope of 
0.5, the equations for which are : 


for = 80-2, Hk, = 0-095 [1 + (x 
for = 180, Hk, 


(11) 
(12) 

The lack of correlation for the various lines is 
clearly due to variations of the iodine concentration. 
The parameter « is a function of the catalyst concen- 
tration C,.and may be a function of liquor rate. When 
graphed in Fig. 7, x is seen to vary linearly with 
C,. The intercept (the value of « when C, = 0) is 


0-164 [1 + (2 


Increase in liquid film coeflicient due to chemical 


increase is proportional to the 0-5 power 
of the same ratio. Equation (3) leads 
to the conclusion that increasing the 
reaction rate constant k, would increase 
the exponent so that it approaches 
unity. However the increase in reaction 
rate due to catalytic action did not alter the value 
of the index. It may be that the index is a function 
of diffusivities. 

Variations in the reaction rate constant have a 
similar effect upon the increase in liquid film coeffi- 
cient as have variations in the concentration of the 
liquid phase reactant. 


® [ =60-2 lb/(hrilft) 
[=180 


08 
107C¢ Lodine concentration 


OF 


Fig. 7. The variation of « as a function of iodine concentration 
Ce. 
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NOTATION 


A = interfacial area, ft.” 
C = concentration of solute, Ib. mol per ft. of solution, 
D = diffusivity, ft.?/hr. 
H = Henry's law coefficient, Ib. mol /(ft.) (atm.). 
K = overall absorption coefficient. 
k = individual film coefficient. 
k’,, the liquid film coefficient for physical absorption. 
k, = reaction rate constant, (hour)- (Ib. mol/ft.4)-" 
Ny = rate of absorption of chlorine, lb. mol/hour. 
P = pressure, atm. 
P py the logarithmic mean of the partial pressure of air 
at the phase boundary and in the bulk of the gas stream. 
AP, _. the logarithmic mean driving force, atm. 
Re = Reynolds number. 
Se = Schmidt number. 
vo = gas velocity, ft./sec. 
2 = the distance of an element of the liquid film from the 
interface, ft. 
a = a parameter defined by equations (11) and (12). 
8 = a constant in equation (9). 
I’ = liquor rate, Ib./(hr.) (ft.). 
» = absolute viscosity, Ib. mass /(ft.) (hr.). 
p = density, Ib./ft.* 
pu = density of the gas stream, Ib. mol per cu. ft. 


Subscripts 
A = refers to chlorine. 
B = refers to 2-ethyl hexene-1. 
C = refers to iodine. 
G = refers to the gas. 
i = refers to the dissolved chlorine at the interface. 
L = refers to the liquor. 
The units used are feet, hours, pounds and atmospheres: 
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The kinetics of the addition of chlorine to olefins in carbon tetrachloride 
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Summary —For the interpretation of the rates of absorption of chlorine from air by solutions of olefins 
and iodine in carbon tetrachloride, it is necessary to know the kinetics of the addition of chlorine to the 
olefins under the conditions of the absorption. The rate of addition of chlorine to 2-ethyl hexene-1 and 
to oleic acid in carbon tetrachloride was found proportional to the concentration of the chlorine and to 
the concentration of the olefin. With excess cyclohexene, the rate of addition of chlorine to cyclohexene 
in carbon tetrachloride was found to be proportional to the second power of the chlorine concentration 
and independent of the concentration of cyclohexene. 


Résumé—Pour interpréter la cinétique du lavage des mélanges chlore-air avee des solutions dans CCl, 
doléfines et d’iode, il est indispensable de connaitre la cinétique de l'addition du chlore aux oléfines dans 
les mémes conditions. 
La vitesse d’addition du chlore a l'ethyl-2 hexéne-1 et & l'acide oléique dissous dans CCl, est propor- 
tionnelle a la concentration du chlore et a celle de l'oléfine. Dans la cas du cyclohéxene, la vitesse d’addition 
- est proportionelle au carré de la concentration du chiore et indépendante du celle du cyclohexéne, si ce 
953 dernier reste toujours en excés par rapport au chlore absorbé. 


INTRODUCTION 


In a previous publication [1] the rates of absorption 
of chlorine from air by solutions of 2-ethyl hexene-1 
in carbon tetrachloride were reported. The rate of 
reaction between chlorine and 2-ethyl hexene-1 in 
solution appreciably affects the rate of absorption, 
and in order to interpret the results of the absorption 
experiments it is necessary to know the kinetics of 
this addition. This paper describes the investigation 
into the kinetics of this addition, as well as the 
addition of chlorine to oleic acid and cyclohexene 
under similar conditions. The rates of absorption of 
chlorine into solutions of these olefins is at present 
being studied. 

The addition of halogens to olefins in non-polar 
solvents has been reviewed by pe La Mare [2] and 
by Remick [3]. Both agree that the additions are 
difficult to reproduce and that the kinetics are often 
intractable. The rate of addition is increased by 
polar molecules -such as hydrogen bromide and 
water —and is decreased by oxygen. The rates of 
addition which are sufficient to influence considerably 
the rates of absorption are much faster than the 
rates of reaction in solution which have been previously 
measured. 

The rates of addition studied in this paper are so 
great that the usual unsteady-state method of 


measurement is unsuitable. A steady-state method 
is required. HovuGen and Watson [4] give a mathe- 
matical treatment of flow reactions for the cases when 
the stream of reactants moves progressively through 
the reactor without longitudinal mixing and without 
appreciable mixing in the direction of flow. These 
conditions are satisfied in a tubular reactor. HouGeN 
and Watson give the reaction volume, V, for a 
specified conversion, 2, for a constant mass feed rate, 
F, as 


where r is the rate of reaction, moles of component 
A converted per unit volume of the reacting system 
per unit time. For a reacting system, the density 
of which is independent of the conversion zx, the 
particular solution for the specific reaction rate, k,, 
for a second-order irreversible reaction is 


F In ©40°R 


k= 
V p (ego — Cao) 


(1) 


CBo 


or, for equal concentrations of the reactants, 


k _ F (€49 — a) 


Vp CAQCA 
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EXPERIMENTAL APPARATUS AND PROCEDURE 
The general arrangement of the apparatus is shown 
in Fig. 1, and the details of the glass tubular reactors. 
The left-hand liquid feed vessel held a solution of an 


contact 


Fig. 1. The general arrangement of the apparatus. 


olefin in carbon tetrachloride. The liquor from cach 
vessel passed downwards through a flowmeter and 


Table 1. The dimensions of the glass tubular reactors 


Reactor No. 1 | 2 | 3 


3 
| 
Length, em 22-9 | 11-4 6-8 13-3 80 
Diameter, cm 0-249 0267 | 0-280 O122  O124 
Volume, ml oh 0-637 | 0-418 0-180 | 0-096 


entered the reactor, a glass tee, from opposite sites. 
The reacting mixture passed downwards from the 
tee through a contact solution (5 normal hydrochloric 
acid) to a product receiving vessel. For each run a 
mass balance was calculated from the analyses of the 
liquor in cach of the liquid feed vessels, the rate of 
feed from each vessel, the analysis of the product, 
and the total liquor rate. If the balance was more 
than 6%, out, the run was rejected. The initial 
concentration of a reactant, ie. ¢4, (chlorine) or 


Co (Olefin), was calculated from its concentration in 
a liquid feed vessel, and the fraction of the total 
liquor rate contributed from this vessel. The con- 
centration of chlorine in the discharge from the 
reactor was calculated from its concentration in the 
contact solution and a distribution coeflicient for 
chlorine between carbon tetrachloride and 5 normal 
hydrochloric acid. The distribution coeflicient was 
determined as a function of temperature by running 
only liquid from the right-hand feed tank (i.e. without 
olefin) and determining the concentration in both 
phases. The value of V used in equations (1) and (2) 
allowed for one-half of the volume of the liquid 
stream as it passed through the contact solution 
assuming that the stream cross-section remained 
constant during this period. Thus all the variables 
appearing in equations (1) and (2) were known. 

Because of the possible effects of water and oxygen, 
the solutions were made from carbon tetrachloride 
saturated with both water and air. The addition 
was carried out in the dark. 


Reaction Products 


Analysis of the reaction products showed that there 
was no hydrogen chloride formed, indicating that no 
substitution was taking place. Within the limits of 
accuracy of the analyses, one mole of chlorine reacted 
with one mole of olefin. 


Materials 


The carbon tetrachloride was purchased from Imperial 
Chemical Industries of Aust. and New Zealand Ltd. of 
Botany, and was labelled “ sulphur-free.” It was 
redistilled in a glass laboratory still, the column of 
which contained twelve theoretical plates, at a reflux 
ratio greater than 10:1. The first and last 10% 
of the distillate was discarded. As collected and used, 
the carbon tetrachloride conformed to the following 
tests of the British Pharmacopocia ; boiling point, 
specific gravity, and the limit tests for ionisable 
chlorides and free chlorine. 

The 2-cthyl hexene-1 was prepared from 2-ethyl 
hexanol and sulphuric acid. The 2-ethyl hexene-1 
as used conformed to the following tests : 

1. Boiling point, 120 121°C, 

2. Specific gravity (15 C/15°C), 0-715 0-717. 

3%. lodine number 225. corresponding to a purity 

of 99-5'Y,. 
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4. When 10 ml were shaken with 10 ml of water, 
and allowed to separate, the aqueous layer was 
neutral to litmus, indicating no free acid or 
alkali. 

The oleic acid was obtained from New South Wales 
Government Stores and conformed to the requirements 
of the British Pharmacopoeia. 

The cyclohexene was prepared from cyclohexanol 
and phosphoric acid. The cyclohexene as used con- 
formed to the following tests : 

1. Boiling point, 82-8-83-6C. 

2. Specific gravity (20°C 4°C), 0-812. 

3. Iodine number 306, corresponding to a purity 

of 99.0%. 

4. No free acid or free alkali. 

The chlorine was supplied by J.CJ.A.N.Z. of 
Botany in steel cylinders. It was not less than 
99.5%, as Cl. It was used as received. 


Resuuts Discussion 
The Addition of Chlorine to 2-Fthyl Hexene-1 


The first reaction chosen for study was between 
chlorine and 2-ethyl hexene-1 in carbon tetrachloride. 
In all the runs, the initial concentration of the olefin 
(Cyc) exceed the initial concentration of the chlorine 
(C49) because this condition corresponds to the 
conditions during the gas absorption experiments. 
The data from the rate measurements are given in 
Table 2. 

The temperatures at which the rates were measured 
fall into three groups ~ nine measurements at 19-4°C 
to 22-8°C, six measurements at 37-8°C to 38-9°C, 
and ten measurements at 12-8°C to 13-9°C. The 
calculated values of the specific reaction rate within 
each group is satisfactorily constant when calculated 
from equation (1). The reaction is thus second-order, 
its kinetics being expressed by 


— (3) 


The data from Table 2 has been grouped according to 
temperature to give the following average results. 


Temperature, °C k, 
21 152 
38 173 
13 132 


The activation energy calculated from these grouped 


results is 2,000 cal. per g. mole. The effect of tempera- 
ture on the specific reaction rate is approximately 
the same as its influence on the diffusivities of the 
reacting molecules. 


The Addition of Chlorine to Oleic acid 


The rate of this addition was measured in the same 
manner as the rate of addition of chlorine to 2-ethyl 
hexene-1. The data for these measurements are given 


Table 2. The addition of chlorine to 2-ethyl hexene-1 


Run | Temp. Vp/F 
No. °C C40 C4 sec. ka 
1 19-4 0-0091 0-0209 0-0023 0-600 152 
2 19-4 0-0091 0-0209 0-0032 0-376 161 
3 19-4 0-0091 0-0209 0-0045 0-284 136 
4 211 0-0110 | 00211 | 00034 | 0525 139 
5 21-1 0-0110 0-0211 0-0046 0-324 180 
6 21-1 00110 | 06-0211 0-0059 0-232 | 157 
7 22-8 0-0160 0-0194 00-0056 | O-GI8 | 136 
8 | 228 | 00160 | 00194 | 00068 0-376 | 168 
9 22-8 0-0160 | 060-0194 | 60-0092 | 0-259 | 143 
10 38-3 0-0165 | 00228 | 0-0045 | 0-586 | 150 
ll 37-8 0-0165 | 0-0228 | 00060 | 0-358 | 180 
12 38-9 0-0165 0-0228 0-0072 0-257 190 
13 38-3 0-0143 0-0242 0-0059 0-241 188 
14 37-8 0-0143 | 00242 | 00052 | 0-339 | 158 
15 38-3 0-0143 0-0242 0-0040 0-548 163 
16 13-9 0-0189 0-0242 0-0059 0-579 128 
17 13-9 0-0189 0-242 0-0079 0-363 141 
18 13-9 0-0189 0-0242 0-0095 0-251 150 
19 13-9 0-0189 | 00242 | 0-0153 | 0-094 | 121 
20 13-9 0-0189 0-0242 0-0170 0-054 | 132 
21 12-8 0-0208 0-0242 0-0181 0-054 115 
22 12:8 0-0208 0-0242 0-0162 0-093 | 114 
23 12-8 0-0186 0-0242 0-0097 0-258 130 
24 12:8 0-0186 0-0242 0-0072 0-350 156 
25 12:8 0-0186 0-0242 0-0056 0-573 135 


in Table 3. The temperatures at which the rates were 
measured fall into two groups - five measurements at 
19-4°C to 23-3°C, and five measurements at 40-0°C. 
to 42-2°C. The calculated values of the specific 
reaction rate in each group is satisfactorily constant ; 
the reaction is therefore second-order, its kinetics 
being expressed by equation (3). The values of the 
specific reaction rates for the addition of chlorine to 
oleic acid are the same, within the limits of error of 
the method, as the values obtained for the specific 
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reaction rate for the addition of chlorine to 2-ethyl 
hexene-1. 


Table 3. The addition of chlorine to oleic acid 


Run | Temp Vp/F 
No. *c © 40 C4 sec. ka 
1 29-6 0-0071 0-0458 0-0030 0-130 | 153 
2 194 0-0076 0-0268 0-0033 O-212 | 162 
3 22-2 0-0205 0-0259 0-0038 0-695 | 170 
4 23-3 0-0124 0-0350 0-0014 0-570 | 140 
5 23-3 0-0186 0-0282 0-0020 0-826 | 169 
6 40-6 0-0091 0-0170 0-0018 0-682 | 195 
7 40-0 0-0064 0-0183 0-0015 0-504 | 191 
8 40-0 0-0064 0-0183 0-0023 0-324 | 200 
9 419 0-0064 0-0183 0-0036 0-169 | 199 
10 42-2 0-0064 0-0183 0-0042 0-103 | 176 


The data from Table 3 has been grouped according 
to temperature to give the following average results : 


Temperature, °C k, 
22 159 
41 192 


The activation energy is the same as that of the 
addition to 2-ethyl hexene-1 within the limits of error 
of the method. 


The Addition of Chlorine to Cyclohexene 


The measurements of the rates of addition were 
carried out under the same conditions. The data is 
given in Table 4. The temperatures at which the rates 
were measured fall into three groups — eight measure- 
ments at 26-7°C to 27-8°C, five measurements at 
18-9°C, to 21-1°C, and five measurements at 35.0°C. 
to 87-8°C. The calculated values of the specific 
reaction rate using equation (1) vary greatly. That 
the reaction product has no influence on the rate of 
addition is shown by run 8, to the feed of which was 
added sufficient 1, 2-dichlorocyclohexene to give an 
initial concentration of 0-130 g. moles per litre. In 
runs 12 and 13, the initial conceatration of the 
cyclohexene was 0.0910 g. moles per litre, a con- 
centration considerably higher than used in the 
other runs. The specific reaction rate, k,, is cor- 
respondingly low. It is possible that the reaction is 
zero order in respect to the cyclohexene. Specific 


reaction rates were calculated, for zero order in 
respect to the olefin, assuming that the rate was 
first, second, or third order in respect to the chlorine 
concentration, The calculated values of the specific 
reaction rate were most nearly constant for the 
assumption of second order in respect to the chlorine 
concentration and zero order in respect to the cyclo- 
hexene concentration. The variation in this specific 
reaction rate, k,, is rather greater than the variations 
previously obtained for k, (i.e. for 2-ethyl hexene and 
oleic acid). The kinetics for the reaction may be 
expressed by 


=keg (4) 


The data from Table 4 has been grouped according 
to temperature to give the following average results : 


Table 4. The addition of chlorine to cyclohexene 


Run | Temp Vp/F 
No. °c C40 sec. | hy 
1 26-7 | 0-0075 | 0-0341 | 0-0025 | 0-547 65 | 500 
2 27-2 | O-O110 | 06-0202 | 0-0032 | 0-203 | 182 | 755 
3 27-8 | 0-0138 | 0-0268 | 06-0063 | 0-102 | 336 | 845 
4 | 26-1 | 60-0051 | 0-0225 | 0-0017 0-568 | 93 | 690 
5 | 267 | 0-0059 | 0-0212 | 0-0028 | 0-248 | 150 | 770 
6 26-7 | 0-0122 | 0-0167 | 0-0055 | 0-161 | 390 | 620 
7 26-7 | O-O118 | O-0170 | 0-0051 | 0-148 | 435 | 750 
8 27-8 | 0-0147 | 0-0185 | 06-0051 | 0-197 | 430 | 650 
9 | 21-1 | 60-0055 | 0-0292 | 0-0014 | 0-875 91 | 610 
10 21-1 | 06-0074 | 0-0201 | 0-0018 | O-721 | 118 | 590 
il 21-1 | 60-0074 | 0-0201 | 0-0017 | 0-839 | 108 | 550 
12 18-9 | O-O112 | 06-0910 | 0-0042 | 0-297 38 | 500 
123 194 | O-O112 | 06-0910 | 0-0013 | 1-49 17 | 480 
14 | 35-0 | 0-0062 | 0-0168 | 0-0034 | 0-147 900 
15 | 35-0 | 0-0067 | 0-0159 | 0-0012 | 0-805 850 
16 36-1 | 0-0069 | 0-0156 | O-0011 | 0-832 920 
17 37-8 | 0-0045 | 0-0202 | 0-0027 | 0-187 790 
18 | 37-2 | 0-0070 | 0-0178 | 06-0016 | 0-670 7 

Temperature, °C hk, 

27 700 

20 550 

36 S40 


From the data, the activation energy is 2,000 cal. per 
g. mole, the figure obtained for the two previous 
reactions. 
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Some independent evidence is available to confirm 
the finding that the rate of addition of chlorine to 
cyclohexene is zero order in regard to the concentra- 
tion of the olefin. The rate of absorption of chlorine 
by solutions of 2-ethyl hexene-1 is increased by 
increasing the concentration of olefin in the liquor. 
Preliminary tests have indicated that oleic acid has 
a similar influence on the rate of absorption of 
chlorine, but that the rate of absorption of chlorine 
is not increased by increasing the concentration of 
cyclohexene in the liquor, provided that there is 
sufficient cyclohexene in the liquor to react with all 
the chlorine absorbed. 
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NOTATION 


¢4 = the concentration of chlorine in the discharge from the 
reactor g moles per litre 
¢ 49 = the concentration of chlorine in the feed to the reactor 
g moles per litre 
Cp = the concentration of olefin in the discharge from the 
reactor g moles per litre 
Cho = the concentration of olefin in the feed to the reactor 
g moles per litre 
F = the total liquor feed rate to the reactor g per sec. 
k, = the specific reaction rate, defined by equations (1), (2) 
and (3) 1./(g mole) (sec.) 
ky = the specific reaction rate, defined by equation (4) 
1. /(g mole) (sec.) 
t = time sec. 
p = density g per ce 
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Overall heat flux values from condensing steam to boiling liquid 
E. A. Jones and L. W. J. Loveess 
Courtaulds Ltd., Research and Development Section, Coventry 
( Received 5 October 1952) 


Summary—From values of the overall heat transfer coeflicients determined over a range of 
temperature differences the corresponding total heat flux values have been calculated for the system 
studied. The temperature difference corresponding with the maximum heat flux is deduced. 


Résumé —L'auteur mesure le coefficient global de transmission thermique pour tout une gamme 


de difference de temperature. 


Il en deduit les valeurs du flux thermique total et la valeur de la 


difference de temperature correspondant au maximum de flux. 


INTRODUCTION 
To check the calculated design heat exchange surface 
of an evaporator calandria and also to obtain some 
experience with the alloy steel tubes proposed, it was 
considered advisable to undertake some small! scale 
tests before proceeding with the manufacture of the 
full scale equipment. 

The proposed evaporator calandria was to be fitted 
with vertical 18/8/3/1 Cr/Ni/Mo/Ti stainless steel 
alloy tubes having an internal diameter of 1-50 inches. 
The small scale apparatus for carrying out the required 
tests was therefore designed to incorporate a length 
of the same type of alloy steel tube. 

Apart from the mechanical performance of the alloy 
steel tube, which was quite satisfactory, the results 
obtained were of particular interest in that they 
clearly showed a maximum value of the total heat 
flux when this was plotted against the mean tem- 
perate difference between the heating and heated 
fluids. 

APPARATUS 


A diagram of the apparatus is given in Fig. 1. 


The evaporator consisted of a vertical length of the 
alloy steel tube I3in. long surrounded by a 3in. 
diameter carbon steel steam jacket of 12in. overall 
length. The evaporator tube was provided with a 
liquid level glass into the top of which the fresh liquor 
feed was metered through a sealing tube. The feed 
liquor was introduced in this manner in order to pre- 
vent crystallisation in the level glass. The top and 
bottom of the evaporator tube were connected with 
similar alloy steel tubes of I}in. and lin. nominal 
diameters respectively to a separating chamber 
which was a simple vessel 18in. high and 6in. 
diameter also made of the same grade of alloy 
steel. 

At the lowest point in the lin. liquor line from 
the separator to the evaporator a drain valve was 
fitted for withdrawal of residue liquors. 

Vapour from the top of the separator passed 
through a lin. diameter copper pipe to a submerged- 
coil type condenser having four coils of lin. diameter 
copper pipe each coil being Tin. in diameter. The 
distillate was collected from the bottom of the 
condenser. 

A 0.5in. diameter thermowell was inserted down 
the centre of the evaporator tube in order to obtain 
the boiling temperature of the mixture. A second 
thermowell was installed in the separator in such a 
way as to indicate the liquid temperature in the 
bottom of this vessel. 

Saturated steam was used for heating the evaporator 
and the pressure of this was measured by a gauge 
fitted as near to the evaporator steam jacket as 
possible. Steam was manually controlled by a needle 
valve. 

The heat losses from the liquor return line from the 
separator to the evaporator were minimised by 
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efficient lagging and by steam tracing with a }in. 
diameter copper pipe carrying 30 psi gauge steam. 


EXPERIMENTAL 


The liquor under investigation consisted of a mixture 
of potassium acetate, acetic acid and water. An 


initial quantity of this liquor was put into the appara- 
tus so that approximately one third of the evaporator 
tube was submerged. Heating was then commenced 
and as soon as proper boiling was initiated the feed 
liquor was introduced at a steady and metered rate. 
Distilled aqueous acetic acid was collected from the 
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Table 1. 


condenser while the total quantity of material in the 
system was maintained constant by controlled with- 
drawal of the residual liquor from the drain valve 
which had been installed for this purpose. 

During the tests the thermometer in the evaporator 
tube gave readings which were somewhat erratic but 
which were generally similar to those indicated by 
the thermometer in the separator. Readings of the 
evaporator liquor temperature were therefore dis- 
continued in favour of the separator liquor tempera- 
ture. 

A number of tests each of approximately three 


Feed rate Ibs. /hour 


Steam pressure psi gauge 


Steam temperature °F. 


Temperature of oo from 
separator °F. 


274 274 


Mean temperature difference °F. 


Distillate rate Ibs. /hour 


Distillate composition wt. 
acetic acid 


Sensible heat supplied to feed 
Btu/hr. 


Latent heat supplied to 
distillate Btu /hr. 


4019 


Total heat transferred Btu/hr. 


Overall heat flux 
Btu/hour/sq. ft. 


Overall heat transfer coefficient 
Btu/hour/sq. ft./°F. 


| 
| 
Trial No. 1 2 3 | 4 5 6 7 8 9 10 
16-2 17-6 16-5 | 5) 17-3 17-4 16-5 16-0 15-1 
Td 68 71 4 91 7 | 96 97 100 111 123 
Ts B15 317 334 332 336 335 336 338 B45 352 
2s0 
an, 39 46 53 35 56 61 64 70 71 72 
392) 418 4 502) 5 473) 4-67) #73 
10-5 10:3 10-7 115 113 11-3 110 115 
Po 2620 | 2790 | 2740 | 2680 | 2870 | 2780 | 2770 | 2580 | 2570 | 2490 
| sss | | | | mmm | sone | | | 
Ps | 5851 | 6234 | 6721 | 6542 | 6641 | 6799 | 6808 | 6362 | 6298 | 6224 
14900 15850 [17100 | 16620 |16900 |17300 |17320 16200 | 16010 
| 220 | 220 
a3 
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hours duration were carried out and for each test 
the feed temperature, temperature of the liquor 
leaving the separator, steam pressure and rates at 
which materials entered and left the system were 
carefully recorded. 


CALCULATION OF RESULTS 


In order to calculate the effective overall heat transfer 
coefficients from the condensing steam to the boiling 
liquid during the periods of steady evaporation, it 
was assumed that the total heat passing through the 
heat transfer area was equal to the sum of the sensible 
heat required to raise the temperature of the feed 
to the temperature of the liquor leaving the separator 
plus the latent heat removed in the distillate, i.e. it 
was assumed that the heat losses from the separator 
and pipes connecting it to the evaporator were 
negligible and that the boiling temperature of the 
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Fig. 2. Overall heat transfer coefficient vs. mean temperature 
difference 


liquid mixture was the same as the temperature of 
the liquid leaving the separator. Both of these 
assumptions affect the quantitative accuracy of the 
results but, since the results from the small scale 
apparatus were unlikely to be identical with those 
obtainable in the full scale equipment, it was consi- 
dered that the assumptions were justifiable and no 
attempt was made to determine the heat loss from 
the apparatus. 

Data employed in calculating the effective overall 
heat transfer coefficients and total heat fluxes were :— 

(a) Area of heat transfer surface taken as 0.393 

sq. ft. 
(6) Latent heat of acetic acid taken as 172 Btu/Ib. 


(c) Average temperature of feed liquor taken as 

60°F. 

The averages of the experimentally observed 
readings together with the values of the overall heat 
transfer coeflicients and heat fluxes calculated from 
them are given in Table 1. The values of the transfer 
coefficients and heat fluxes are plotted against the 
mean temperature differences in the Graphs Figs. 2 & 3. 
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Fig. 3. Heat flux vs. mean temperature difference 


Discussion oF THE RESULTS 


Figs. 2 and 3 show that the maximum heat flux 
value obtainable in the apparatus in which the tests 
were made was 17,350 Btu /hr/sq. ft. and that this 
occurred when the mean temperature difference was 
approximately 62°F. The overall heat transfer 
coeflicient corresponding with this maximum value 
of the heat flux was 275 Btu hr/sq. ft./°F. 

Based on the results of these tests the most econo- 
mical design of evaporator calandria for the particular 
service under consideration would clearly be the one 
utilising a mean temperature difference of approxi- 
mately 62°F. 

The graphs also demonstrate clearly that an increase 
in steam pressure, i.e. an increase in the mean tem- 
perature difference, is not necessarily synonymous 
with an increased rate of boiling and may in fact 
actually operate in the reverse direction when the 
mean temperature difference corresponding with 
maximum heat flux has been exceeded. 
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A comparative study on the rate of mixing in stirred tanks 


H. Kramers, G. M. Baars and W. H. Kwnou 
Laboratorium voor Physische Technologie, Technical University, Delft, Netherlands 
(Received 6 November 1952) 
Summary— Measurements have been carried on the rapidity of distribution of a small amount of liquid 


added to the contents of an agitated vessel. 


With the adopted measuring method mixing speeds under 


different conditions can be compared. Results are given for two different tank sizes, three kinds of 
agitators at different locations and speeds, with and without baffles. 
Résumé—Des mesures ont été effectuées pour déterminer la vitesse de distribution d'une faible quantité 


de liquide qui est ajoutée au contenu d'un réservoir muni d'un agitateur en mouvement. 


La méthode 


de mesure appliquée permet la comparaison des vitesses de mélange sous des conditions variées. Les 
résultats donnés portent sur deux dimensions différentes du réservoir et trois espéces d’agitateurs qui sont 
disposés & des endroits différents ct opérent & des vitesses variées, avec ou sans chicanes. 


1. INTRODUCTION 


Few chemical engineering operations are as evasive 
to experimental investigation as those connected with 
mechanical agitation. The choice and the design of 
equipment to be used for one of the many specific 
objectives for which agitation is required depends 
largely on experience. Only the power consumption 
of a number of standard types of agitators can be 
predicted within a reasonable margin thanks mainly 
to the work published by Ruswrown et al. [1]. 

This paper deals with batch mixing of two liquids 
in a cylindrical vessel which is provided with a 
stirring device. A fundamental approach to this 
particular problem would imply the measurement of 
the average velocity of the stirred liquid and the 
degree of turbulence at any place in the vessel, 
because essentially the problem is a hydrodynamic 
one. But in view of the experimental difliculties 
involved it is only natural that an investigation into 
the effect of such a flow pattern seems to be more 
attractive, although the results may contribute less 
to a general understanding of the problem. 

Thus, in the present investigation the completion 
of the mixing process was observed by means of 
concentration measurements. Clearly, a complete 
record of the distribution of the concentration 
throughout the vessel as a function of time would 
give the most complete information on the mixing 
performance of the system. This again is prohibitive 
for practical reasons, so that we contented ourselves 
to measuring the concentration in two fixed points, 
appropriately chosen, and deriving from these an 
arbitrary quantity, the mixing time ©, which is 
needed to obtain a certain degree of uniformity of 


the mixture. The results thus obtained are of com- 
parative interest only and no general relationships for 
design purposes can be expected from them. In 
mixing, however, even comparative data are so 
searce that results of this nature may be of sufficient 
interest. 


2. PrincipLeE oF MEASUREMENT OF THE MIXING 
Time 


The mixing vessel, which contains a weak electrolytic 
solution, is operating at the circumstances under 
investigation. At a certain time a small amount of 
a concentrated solution is thrown into the liquid. The 
effect of this disturbance is measured by means of 
two measuring cells and a record is made of the 
difference in concentration between those two. As 
the mixing proceeds the average value of this differ- 
ence gradually decreases and the mixing time @ is 
defined as the time interval between the injection of 
the disturbance and the moment from which on the 
measured concentration difference remains less than 
0-1 per cent of the average concentration in the 
vessel, 

Woop, Wuirreman and Bapcer [2] and also 
Grese [3] used a somewhat similar method in their 
investigations into the mixing effect of paddles. A 
great advantage of the injection method is that the 
stationary flow pattern is practically not disturbed 
by the measurement. 


3. 
The mixing experiments were carried out in two 
cylindrical vessels with diameters of 0.32 and 0-64 m 
respectively. In all cases the liquid height was equal 
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(a) (b) 


to this diameter. For most of the experiments we 
used a 3-blade marine propeller having a diameter 
1/4 of the tank diameter. The two propellers for the 
two vessels were exactly similar, their dimensions 


i 

? 1 

| 
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Fig. 2. Kinds of baffles investigated (dimensions expressed 
as a fraction of the tank diameter D). 


differing by a factor 2. Further particulars are given 
in Fig. 1. This figure also shows the two impeller-type 
agitators, which were used in the large tank only. The 
position of the agitator and the direction of the shaft 


(ce) 
Fig. 1. Agitators used in the investigation (dimensions expressedas a fraction of the diameter L). 


were variable. Its rotational speed could be varied 
from 300 to 2,000 r.p.m. For the large tank means 
were available for measuring the torque on the rotor 
shaft. 

Either no baffles were used or one of the arrange- 
ments as shown in Fig. 2. Those of Fig. 2b have been 
mentioned by Newirrt [4] to produce more turbulence 
than the more usual arrangement of Fig. 2a. The 
crossed baffles on the tank bottom have been suggested 
by Reave.t [5] to be used with a propeller for the 
suspension of particles. 

The injection of the strong solution was performed 
by mechanical means. A measured amount of solution 
was discharged from a small reservoir within a few 
tenths of a second by operating on electromagnetic 
shutter device. The place of injection always was 
half-way between the wall and the axis of the vessel, 
at the liquid surface. We used KCI as the dissolved 
electrolyte because the conductivity of its solution 
in water is proportional to the concentration from 
1 to 20 mg cm*. 

The two similar cells for measuring the electric 
conductivity consisted of a few turns of Pt-wire, 
forming a cylinder of 25 x 25mm, as the ground 
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electrode, and a Pt-wire in the axis of this cylinder 
as the other electrode. Both cclls formed part of a 
Wheatstone resistance bridge in which additional 
compensation for capacity of the cells was possible. 
The 1,000 c.p.s. cross-potential was amplified, passed 
through a phase sensitive rectifier and fed to a fast 
recorder (Brush Development Co.). To protect the 
recorder against overload, due to relatively high 
output voltages occurring immediately after injection 
of the strong solution, the maximum amplitude of 
these voltages was limited by the special rectifier 
circuit. Fig. 3 shows a block diagram of the measuring 


7 volts Lt 


1000 cps 


Amplifier 


Rech fier 


‘onductivil | | 


Mixing tank Recorder 


Fig. 3. Block diagram of measuring arrangement. 


system. In Fig. 4 a few records are shown of the 
difference in electrical conductivity between the 
two measuring points as a function of time during 
a mixing experiment. Apart from this, the moment of 
injection was recorded and after completion of the 
mixing a constant deviation of 0-1 per cent of the 
average concentration was marked. This was achieved 
by short-circuiting a small part of one of the fixed 
resistances (R,) in the bridge. Fig. 4 clearly indicates 
how the “mixing time” 6 was read from the record 
taken. 

The criterion for sufficient mixing at 0-1 per cent 
relative deviation is of course arbitrary, as well as 
the other measuring conditions. In order to get 
comparative results, the latter were standardized 
according to the following specifications : 

Liquid in tank: tap water with 2.5 to 3 mg/cm® 

KCl, room temperature. 


Injection liquid: V’,, = 10cm? (small tank) or 
80 cm® (large tank) with a concentration ¢, = 
100 mg/em® KCI. 

Size of cells: as stated above, for both tanks 
because for nearly sufficient mixing the dimen- 
sions of the fluctuations are large compared with 
the cell dimensions. 

Position of cells: 1/8th of tank diameter above 
the bottom and under the liquid surface. 

These specifications were arrived at by a great 

number of preliminary experiments ; their significance 
will be treated in Section 5. 

From the records in Fig. 4 it can be seen that a 

considerable spread in the values of @ was to be 


HH: 


Fig. 4. Two mixing records obtained under the same circum- 
stances, indicating how the mixing time 6 was determined. 


expected, due to the fact that on a signal which is 
gradually decreasing to zero random fluctuations of 
the same order of magnitude are superimposed. From 
about 700 readings we found the distribution of @ 
to be fairly normal with a standard deviation of about 
15 per cent. We decided to repeat each measurement 
15 times. Thus the average 6 values, as used for the 
results lie with 90 per cent probability within + 6-4 
per cent from the actual mean value. A difference of 
more than 10 per cent between two @ values can be 


considered to be highly significant. 


4. EXPERIMENTAL RESULTS 
The influence of a number of variables on the arbi- 
trarily defined mixing time 6 were investigated : 
(a) rotational speed of the agitator (n) ; 
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(6) direction of rotation, for propeller only ; 
(c) location of the propeller and inclination of the 
shaft ; 
(d) presence of baffles ; 
(e) type of rotor. 
(a) From 8 special runs it was found that @ was 
inversely proportional to the rotational speed of the 


agitator. Fig. 5 gives a few examples of the results. 
20 
sec 
1s Z 
C= 
© without baffles 
Wa | 
2007000800 600 S00 
rpm 
0 005 07 


Fig. 5. Examples of proportionality between 6 and 1/n ; pro- 
peller with vertical shaft, small tank. 


In all these cases the product n @ appeared to be 
constant well within + 10 per cent. Subsequently 
all measurements under conditions (b) to (e) inclusive 


were taken at only two values of n, as much different 
as possible. A statistical analysis of these measure- 
ments confirmed the conclusion given above. Thus 
it was possible to express the efliciency of mixing 
in terms of the product n 6. For this we took the 
number of revolutions N for “sufficient mixing” under 
the standard conditions mentioned above. Since @ is 
expressed in sec and n in rev/min, we have 


N =n 0/60. 


All experiments were carried out in the flow region 
45,000 < Re < 200,000, where Reynolds number is 
defined in the usual way : 


So the constancy of N for a given geometrical 
configuration has only been established between those 
limits. It is likely, however, that also for higher 
Reynolds numbers N will remain constant, because 
for these high speed agitators turbulence has nearly 
fully been developed in the neighbourhood of Re =10*. 
This can be deduced from the measurements of Rusu- 
TON et al. [1] who showed that from that Reynolds’ 
value on the power number N, changes only little, 
at least for turbines and propellers in baffled tanks. 
(+) The measurements in the small tank were taken 
for both directions of rotation of the propeller which 
was always located in the lower half of the liquid. 
Neither appreciable, nor systematic differences be- 
tween the two cases were found. From the obtained 


Table 1 
Number of revolutions N for sufficient mixing. Marine propellerL = D/4, vertical shaft, no baffles, both tanks. 
db | 3D/8 D/2 Tank size 

0 133 121 130 small 
177 176") large 
| D/8 182 121 117 small 
| 218 161 142°?) large 
D/s 157 135 92 small 

h 145 140 110) 
| | large 
0 5D/16 112 92 90 small 
| 

“)Power number Np = 0-26. Np = 0-34. Np = 0-36. 
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Table 2 
Number of revolutions N for sufficient mixing. Marine propeller L = D/4, inclined position, h = D/2, ¢ = D/4, no baffles, 
both tanks. 
a e N tank size 
30° 0 80 small 
30° 0 89 large 
30° D/4 110 large 
15° 0 large 
15° D/4 110 large 


data, for which we refer to [6], we concluded that in 
general a downward thrust of the propeller was to 
be preferred. The data given below only apply to 
the latter case. 

(c) Table 1 presents the final results for the propeller 
with a vertical shaft in a vessel without baffles. The 
main object of excentric location is the improvement 
of the pumping action of the propeller by reducing 
the intensity of the tangential swirl of the liquid. The 
results show that a sufficient excentricity considerably 
improves the mixing performance provided the pro- 
peller is not too near the bottom of the tank. A small 
excentricity, however, and a short distance between 
the propeller and the bottom appear to be unfavoura- 
ble (e.g. e = D/8 and h = D/4), even compared with 
the case of a central shaft. 


Table 3 
Number of revolutions N for sufficient mixing. 
Marine propeller on central vertical shaft, L = D/4, various 
kinds of baffles, both tanks. 


h = 
baffles tank size 
Pig. D/4 3D/8 D/2 
2a 93 90 91 small 
2a = 100 92 large 
2b — 119 103 large 
2e 107 175 140 small 


Only few measurements have been made with an 
inclined shaft. Table 2 gives the results which show 
no great variations. 

(d) The influence of baffles on the performance 
of a central propeller (Table 3) is similar to that of 
excentric location. It appears to be more favourable 
to have vertical baffles against the wall (Fig. 2a) than 
some distance from the wall (Fig. 2b). The crossed 
baffles on the bottom (Fig. 2c) are only effective for 
mixing if the propeller is placed immediately above it. 

(e) Table 4 shows two interesting features. The 
first is that the agitator models Fig. 2b and Fig. 2c 
can give as rapid mixing as a propeller, but that the 
power consumption under the same conditions is of 


Table 4 
Number of revolutions N for sufficient mixing and power 


number Np. Influence of difference rotor types on central 
vertical shaft, L = D/4, without and with baffles, large tank. 
A= 
rotor | baffles Np for 
Fig. Fig. D/4 | 38D/8 D/2 | h=D/2 
la — _ 177 176 0-26 
la 2a _— 100 92 0-40 
la 2b 119 103 — 
1b 2a 113 118 140 2-7 
Ib 2b 75 80 86 — 
le 2b 63 _- 65 47 
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the order of ten times greater. In the second place, 
with these turbine agitators it is profitable to shift 
the vertical baffles towards the central axis. We did 
not try to find the optimum baffling conditions for 
these agitators. 

The limitations of the measuring method and the 
variety of the experimental conditions make a 
general discussion of these results impossible. We 
will add only a few remarks on 

(f) vortex formation ; 

(g) power consumption ; 

(hk) size enlargement. 

(f) A central propeller mixer in a non-baffled 
tank works under unfavourable conditions because 
of the tangential swirl of the liquid which is only 
checked to some extent by the tank wall. At a 
certain speed of revolution (n,,,) of the propeller it 
starts drawing air and its pumping efficiency sharply 
diminishes. 

For the small tank without baffles n,.. was 
measured as a function of the propeller excentricity 
(Table 5). 

Table 5 


Approximate maximum rotor speed n,,,,. propeller with vertical 
shaft, h = D/4, small tank, no baffles. 


D/s 5D/16 


wall 


It is to be noted that a small excentricity (¢= _D/8) 
increases the tendency towards vortex formation. It 
also involves a greater value of N (compare Table 1). 

From similar data, which will not be reported here 
because of their qualitative nature, we obtained the 
impression that as a general rule a flow pattern which 
shows vortex formation at the liquid surface corres- 
ponds with a relatively low mixing rate. This state- 
ment cannot be reversed, however. 

(g) For several combinations in the large tank the 
torque exerted on the agitator shaft was measured 
with a dynamometer. From this the power consump- 
tion P of the agitator could be calculated, from which 
the dimensionless power number Np was derived : 


P 


Ne = 


In the range of experimental conditions, 40,000 < 
Re < 100,000, N, was found to be independent of 
Reynolds number. The N, values observed have 
been added to Tables 1 and 4. As far as a comparison 
with Rusuton’s data [1] permits, they agree well 
with these. The relatively heavy power demand of 
the turbine models according to Figs. 1b and le for 
obtaining about the same mixing performance as 
obtained with a propeller has already been mentioned 
in (e). 

(h) The two sizes of equipment were chosen in 
order to see whether mixing times for large scale 
equipment could be predicted from smaller scale 
experiments. Comparison of geometrically similar 
cases (Tables 1, 2 and 3) shows that N remains of 
the same order of magnitude as long as the tendency 
to form a vortex is small (by excentricity or baffles). 
Restricting ourselves to this condition we may 
calculate the total energy required for sufficient 
mixing : 


P @ = N,N (n/60)?. 


Since both N, and N are independent of Re for 
great Reynolds numbers and in the absence of vortices 
gravity forces do not play a part, we are led to the 
conclusion that in scaling up we have : 


P = constant pL‘ n?. 
Or for the mixing energy per unit volume : 


This means that for geometrical similarity, where 
the linear velocity of the agitator tips and the liquid 
density are kept constant, also the total mixing 
energy per unit volume would remain constant. Under 
these conditions the mixing time @ would be pro- 
portional to the linear dimension L. 

This general and somewhat tentative conclusion is 
not valid for those cases where vortex formation in 
the liquid surface is observed. From Table 1 it can 
be seen that for a central and slightly excentric 
propeller in a tank without baffles N appreciably 
increases with increasing diameter. Apparently not 
only Reynolds number must be regarded to be charac- 
teristic of the flow conditions, but also Froude number 
(Fr = gL (n,;60)), which quantity accounts for the 
influence of gravitationa! forces. It is to be noted 
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that for this case also Np, depends on both Re and 
Fr [1]. The experimental results are too scarce 
however to warrant any general conclusions for this 
flow condition. 


5. SIGNIFICANCE OF THE MEASURED 
“Mixinc Tike” 6 
As @ has only been arbitrarily defined it is of interest 
to know whether the general results depend in some 
way on the chosen conditions of measurement. Apart 
from the operating conditions the following factors 
are of influence on the measured mixing time : 

(h) size and location of the measuring cells ; 

(7) salt concentration of the liquid in the tank ; 

(j) salt concentration and amount of injected 

liquid ; 

(k) criterion for suflicient mixing. 

(h) No experiments were performed with different 
sizes of measuring cells. With the adopted size not 
only long range concentration differences but also 
fluctuations of the concentrations were recorded. The 
great number of repetitions of a measurement can be 
regarded as an averaging process of the latter fluctua- 
tions. 

As to the location of the cells, the applied measuring 
method has only significance if there is no correlation 
between the concentrations in the two places of 
measurement. Further, from observations with one 
cell only it was found that completion of the mixing 
is first reached in the immediate surroundings of the 
agitator. As the mixing proceeds this homogeneous 
region expands and finally contains the whole con- 
tents of the vessel. So in order not to find too short 
mixing times the cells had to be placed as far apart 
as practically possible. 

(i) The measurements of resistances with a 
Wheatstone bridge involve that concentration 
differences can only be measured relative to the 
average concentration in the vessel. All other con- 
ditions remaining constant, one would find a longer 
mixing time for a low tank concentration than for a 
high one. This was verified in a number of experi- 
ments, on account of which the most convenient 
concentration interval of the tank liquid was chosen. 

(j) The mixing times reported were found by 
introducing a standard disturbance (c,, and V,,) at 
the liquid surface. In order to see whether the same 
relative results could be obtained for other conditions 


of injection, for a few mixing conditions both the 
concentration (¢;) and the volume (V,) of the injected 
liquid were varied. Table 6 shows the experimental 
results obtained in the small tank. They have been 
plotted on a semi-logarithmic scale in Fig. 6 which 
also contains a few similar data for the large tank. 


Table 6 


Number of revolutions N for sufficient mixing. Variation of 
concentration (¢;) and amount (1';) of injected liquid, other 
measuring conditions being according to standards. 


Mixing conditions 

(mg/em*) | (em?) Cie" in 1 2 

5 69 47 
20 10 o2 | 6 42 52 
80 3 O-4 | 104 36 68 
rT) 10 O-4 104 57 71 
20 20 52 70 
80 10 O-8 119 67 a7 
20 14 67 
200 5 10 - 77 3 
100) 10 121 79 
50 20 10 — 73 90 
80 20 16 160 83 oF 
100 20 20 163 85 100 
125 20 2-5 171 92 106 
200 20 40 190 

Q) Cis (2) — Vis (3) see Fig. 6 


It appears that c; and V; do not independently 
influence the mixing time, but only their product 
which is equal to the amount of salt added to the 
contents of the vessel. This is what one would 
expect provided c; is high with respect to the average 
concentration in the tank and V; small compared 
with the volume of the tank. Further Fig. 6 shows 
that N is proportional to log (c,V,), the straight 
lines approximately intersecting near N = 0. This 
means that a different choice of standard disturbance, 
at least within the region investigated, would not 
have altered our comparative results as to the mixing 
performance under different conditions of agitation. 
Finally it can be seen that the slope of the straight 
lines in Fig. 6 are indicative of the rapidity of mixing 
as well as the values of N at one standard condition 
of measurement. However, from this plot it cannot 
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be concluded that the degree of mixing increases 
exponentially with time. 

(k) If one changes the criterion for sufficient 
mixing (in all experiments 0-1 per cent deviation in 
concentration) other mixing times are found. The 
obtained records did not permit a statistical investi- 
gation into the relation between the amount of 
deviation and the time of mixing. From a number 
of separate runs it was ascertained that the same 
mixing time was found for a standard injection and 
a 0-1 per cent criterion as for the injection of twice 
the amount of salt with a 0-2 per cent criterion. 


100 


A-e 


1$0 200 
Fig. 6. Influence of the amount of injected salt on N. 

1. Central propeller in small tank without baffles, A = 3D/8. 

2. Excentric propeller in small tank, ¢ = D/4,h = D/2. 


3. Central propeller in small (vy) and large ( Ww) tank with 
baffles (Fig. 2a), A = 3D/8. 


4. Turbine (Fig. 1c) in large tank with baffles (Fig. 2c), A — D/4. 


In this section we have put forward the limitations 
inherent to the measuring method, at the same time 


showing that the results described in the preceding 
section have a more general meaning provided they 
are considered as relative and not as absolute data. 
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NoTaTION 


= salt concentration of injected solution (mg/cm). 

= tank diameter (m). 

= excentricity of agitator shaft (m). 

= gravitational acceleration (m /sec?). 

= height of agitator above tank bottom (m). 

= diameter of agitator (m). 
rotational speed of agitator (rev /min.). 

/60, number of revolutions for sufficient mixing. 
power consumption of agitator (Nm /sec.). 

; = volume of injected solution (em). 

6 = time for sufficient mixing (sec.). 

vy = kinematic viscosity (m?/sec.). 

p = liquid density (kg/m). 


"S422 > en 


Np = Power number for agitator = 


Fr = Froude number = 


p L3(n/60) 
nL? 
Re = Reynolds number for agitator = =~ 
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Low Temperature Physics —Four Lectures. By F.E.Simon, 
N. Kurti, J. F. and K, Menpeissonn. 132 pp. 58 


illustrations. Pergamon Press Ltd., London 1952. 21s. 


These lectures constituted a course given at the Royal Institu- 
tion in February and March of 1950. Before publication they 
were revised to include some of the more recent developments, 
even up to 1952. 

The titles of the lectures are : 

Low temperature problems ; a general survey (Sion) ; 

The temperature range below 1° absolute (Kurt) ; 

Liquid helium (ALLEN) ; 

Superconductivity (MENDELSSOHN). 

The subject dealt with is thus the phenomena occurring 
at extreme low temperatures (not exceeding 10° K). For the 
chemical engineer this field will only be of theoretical interest, 
since there are (as yet?) no industrial applications. 

However, all those who are interested will very much enjoy 
the brief, simply and clearly written, and well-founded survey 
given of this field, which is of the utmost importance for pure 
physics. P. J. Harincuuizen 


J. TayLor : Detonation in condensed explosives. 196 pp. 
+ 10 plates. Oxford, Clarendon Press 1952. 25s. 


L’auteur traite en 200 pages de tout ce qu'un technicien peut 
désirer savoir sur le phénoméne de la détonation, envisagé 
comme un processus physico-chimique, la question des effets 
des détonations et de leurs applications étant en dehors du 
champ de louvrage. Dans les cing premiers chapitres, qui 
servent d'introduction, auteur présente, avec divers exemples 
et tables numériques & l'appui, des méthodes de calcul des 
caractéristiques des gaz engendrés par l'explosion d'un systéme 
chimique explosif donné, c’est-a-dire la composition des pro- 
duits résultant de lexplosion, leur température, la pressiou 
sous laquelle ils se trouvent ; ces calculs doivent se faire en 
tenant compte du caractére tres imparfait des gaz réels ce qui 
implique d’avoir & passer par les fugacités : toute la difficulté 
est de disposer d'une équation d'état convenable pour ce calcul 
des fugacités ; lauteur a fait choix de valeurs des paramétres 
de cette équation, que l'on pourrait peut-étre discuter mais 
qui conduisent a des résultats satisfaisants. 

L’auteur expose ensuite la théorie classique de CuapmMan- 
Joucuet de la détonation, et-son application d’abord aux 
explosifs ne donnant que des gaz puis a ceux, et c'est le cas 
des explosifs pour mines grisouteuses, dont I'explosion donne 
des produits solides. Des mesures expérimentales de vitesse 
de détonation confrontées aux calculs numériques montrent 
que, dans de nombreux cas, la détonation d'une cartouche 
cylindrique d'un explosif a lieu assez sensiblement en accord 
avec le schéma classique, que l'on peut appeler schéma des 
détonations idéales. Mais il est certains phénoménes, qui 
jouent d'ailleurs un réle important dans diverses applications 
pratiques, qui ne peuvent se comprendre qu’en faisant appel 


a des schémas plus compliqués ; ce sont la dépendance de la 
vitesse de détonation avec le diamétre et l'effet de confinement 
(nature de l'enveloppe de la cartouche): J. TayLtor expose 
les deux théories qui ont été édifiées, au cours de ces dix derni¢res 
années pour en rendre compté. Un chapitre est consacré aux 
deux régimes de détonation (a vitesse élevée et a vitesse faible) 
des >xplosifs liquides ou gélatineux. Enfin un dernier chapitre 
traite de la nature de la zone de réaction, en particulier de sa 
fongueur et de la durée de réaction ; a cet égard, l'étude des 
phénoménes de détonation parait capable d’apporter des 
données précieuses au physico-chimiste, dans le domaine de 
la cinétique des réactions non homogénes a haute température. 
Louts Meparp 


Physikalisches Worterbuch. (Two parts in one volume), 
(Editor: H. Westruar). Publishers: Springer Verlag. 
Berlin, 1952. 833 and 795 pages. DM 148 or £12 19s. 1d. 


The publication of a Dictionary of Physics is an event of 
some importance. No such Dictionary has appeared in English 
since Glazebrook and none in German since the Physikalisches 
Handworterbuch nearly twenty years ago. The present Dic- 
tionary is indeed to be regarded as a successor to the latter 
work, but is intended to cover a wider field, including physical 
chemistry, astrophysics, and to some extent geophysics and 
biophysics. It does not claim to cover applied physics except 
in so far as it is of value for pure research, being intended for 
the physicist who wishes to obtain information on branches of 
physics in which he is not a specialist or for the non-physicist 
who, possessing a certain amount of physical knowledge, needs 
information on the ideas and problems of physics. 

It must be said at once that this book will prove very 
useful : but there are unfortunately several respects in which 
it falls short of the standard which it might have been expected 
to attain. This is particularly true of the treatment of certain 
aspects of modern physics, where some of the articles are clearly 
* second-hand ” accounts. This is perhaps not surprising when 
one remembers the losses that German physics has suffered 
since 1933 and when one sees that of the 80 contributors to the 
Dictionary all but 6 are resident in Germany to-day. To take 
the subject of atomic energy, for example, there is a fairly 
good general description of a nuclear pile but no reference to 
or description of existing piles. There is a very inadequate 
account of neutron spectroscopy with no indication that any- 
thing can be learned about the structure of matter by neutron 
diffraction, let alone any reference to the results obtained. 
There is, too, practically no mention of an atomic bomb, cer- 
tainly not under that heading. Moving farther afield, the 
reviewer could find no mention of antiferromagnetism, multiple 
beam interferometry, interference microscopy or epitaxy, to 
take a few subjects chosen at random; and other subjects 
were not always found to be as complete or up-to-date as they 
might have been. For example, in an otherwise very good 
series of articles on valency no account is given of Pauiinc's 
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recent derivations of valencies from magnetic cousiderations 
although his 1949 Royal Society paper, which deals with his 
latest ideas, is actually included in the bibliography. The 
bibliography itself is such as to suggest that the Dictionary is 
intended only for German-speaking readers, as the works 
quoted are, wherever possible, those published in that language. 
In the article on the Compton effect, for example, there is no 
reference in the bibliography to Compron and AL1ison’s 
standard work on x-rays; and in the majority of instances a 
translation into German or a review account in that language 
is preferred to an original account in another language. No 
bibliography at all is given in the article on biophysics, although 
this deals, among other things, with the effects of radiation on 
living matter, a subject on which there is no lack of published 
work. 
Enough has, perhaps, been said to indicate the way in 
which, in many instances, the treatment of various subjects 
falls short of the desired standard ; and it may be worth while 
turning our attention to the editorial aspects of the work. The 
editing of such a volume as that under review, a volume con- 
taining definitions of about 10,500 terms with about 1,800 
diagrams, and running to over 1,600 pages, is a formidable 
task; and some slips or oversights are inevitable, especially 
in a first edition. There are, however, certain points which 
must be mentioned. First, all the articles are unsigned, which 
is a pity in the reviewer's opinion. Further, in some instances 
cross references are given to terms which are not given as 
quoted. For example, under “ Absorptionsquerschnitt ” 
reference is made to “ Einfangreaktionen,” but this term is 
not listed, and the term actually defined is * Kinfangprozess.” 
Again, the same physical concept is often defined in two or 
more places, and sometimes differently. For example, in the 
definition given under the heading “ Bremsstrahlung” the 
term is restricted to the production of continuous x-rays in 
an x-ray tube, and this point of view is maintained in the article 
on x-rays ; but in that on y-rays the term is rightly extended 
to cover the production of x-rays by a sudden deceleration of 
electrons, although the author of the article seems inclined to 
regard the resulting radiation not as x-rays but as y-rays. 
Another example is provided by the treatment of the electronic 
configurations in the heavy elements. Under the term “Actini- 
den " the modern view is given that the series of elements of 
atomic number 90 and upwards is to be considered as a “ 5-f" 
series, the elements in this series differing in the number of 
electrons in their 5-f shells ; but in Table 11, “ Schalenaufbau 
der Elemente,” the differences appear in the 6-d shell, although 
the configurations are admitted to be uncertain. In this table, 
too, the series stops at uranium. On the other hand, a good 
presentation is given in Fig. 6, under “ Periodisches System der 
Elemente.” No doubt such overlapping and discrepancies are 
very difficult to avoid but they undoubtedly detract from the 
value of the work. 


The Dictionary is arranged strictly alphabetically, which 
appears to be satisfactory. There are several appendices, con- 
sisting of conversion tables of physical units, the Periodic 
table, ete., tables of physical constants, a supplement to the 
dictionary proper, a short history of physics and a list of 
well-known physicists with their dates : the latter list admittedly 
contains mainly the names of physicists from German-speaking 
regions. Not a few misprints occur, but no serious ones were 
noted. 

In spite of the faults of the Dictionary, and they are many, 
it can still be said that it is a valuable book. The reviewer, for 
one, is glad to have it on his shelves. J. THew.is 


Bruno Riepicer; Berechnung von Fraktionierkolonnen 
fur Vielstoffgemische. 50 pp. text, 36 pp. tables. Springer- 
Verlag Berlin, Gottingen, Heidelberg 1951. Price D.M. 18 
or 31s, 6d. 


It is disappointing to see that, in a recent text on fractiona- 
tion, Raouvt'’s law and the thermal properties of pure sub- 
stances are still used to calculate vapour-liquid equilbria and 
the heat of vaporization in a mixture. The author is appa- 
rently not familiar with the enormous amount of research on 
equilibrium constants in particular of hydrocarbons and the 
thermal properties of components when they occur in a mix- 
ture. In Table 15 the vapour pressures of paraffin hydro- 
carbons, defined by their boiling points at 1 kg/em? abs, are 
represented as a function of temperature. Beyond the critical 
temperature of the pure component its vapour pressure has been 
extrapolated by the aid of a special nomogramme analogous 
to the Cox chart, a procedure that one would certainly not 
expect in a modern book on distillation. It is only on the last 
page of the text (p. 50) that fugacities are mentioned, when 
it is said that the author has not investigated whether the 
use of fugacities or activities instead of pressures would be 
of advantage since it is not possible to represent by a simple 
graphic method deviations from ideal behaviour in multi- 
component mixtures ! 

The subjects treated are pressure-temperature curves, 
equilibria of multicomponent mixtures (dew and bubble point 
calculations), flash vaporization, enthalpy of vapour-liquid 
mixtures [the latent heat of the pure substance is used for this 
calculation and it is stated (p. 21) that there are no data about 
latent heat of a component in a mixture at temperatures 
higher than the critical temperature of the pure component!), 
calculations of distillation of binary mixtures according to 
McCase’s and Tutece’s method, and of multicomponent 
mixtures by plate to plate computation. 

The book is certainly not at the level of Hauspranp’s, 
THORMANN's and Kirscupaum’s famous monographs, but it 
may be used as a first introduction. The author's enthusiasm 
for his subject is reflected in the style in which the treatise is 
written. W. R. van Wux 
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